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ABSTRACT 

We report the first detailed chemical abundances for 5 globular clusters (GCs) in M31 from high- 
resolution (i? ^ 25,000) spectroscopy of their integrated light. These GCs are the first in a larger set of 
clusters observed as part of an ongoing project to study the formation history of M31 and its globular 
cluster population. The data presented here were obtained with the HIRES echelle spectrograph on 
the Keck I telescope, and are analyzed using a new integrated light spectra analysis method that we 
have developed. In these clusters, we measure abundances for Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, 
Co, Ni, Y, and Ba, ages >10 Gyrs, and a range in [Fe/H] of —0.9 to —2.2. As is typical of Milky 
Way GCs, we find these M31 GCs to be enhanced in the a-elements Ca, Si, and Ti relative to Fe. 
We also find [Mg/Fe] to be low relative to other [a/Fe], and [Al/Fe] to be enhanced in the integrated 
light abundances. These results imply that abundances of Mg, Al (and likely O, Na) recovered from 
integrated light do display the inter- and intra-cluster abundance variations seen in individual Milky 
Way GC stars, and that special care should be taken in the future in interpreting low or high resolution 
integrated light abundances of globular clusters that are based on Mg-dominated absorption features. 
Fe-peak and the neutron-capture elements Ba and Y also follow Milky Way abundance trends. We 
also present high-precision velocity dispersion measurements for all 5 M31 GCs, as well as independent 
constraints on the reddening toward the clusters from our analysis. 

Subject headings: galaxies: individual (M31) — galaxies: star clusters — galaxies: abundances — 
globular clusters: general 



1. INTRODUCTION 

Stellar atmospheres largely retain the chemical compo- 
sition of the gas from which they formed, and thus con- 
tain a record of the gas chemistry of a galaxy throughout 
its star formation history. The abundances relative to Fe 
of some key elements, a-elements in particular, can be 
used to identify the timescales and rates of star forma- 
tion over the lifetime of a galaxy, a— elements (C, O, 
Ne, Mg, Si, S, Ar, Ca, and Ti) are produced prima,rily in 
type II supernovae (SNII) (e.g. lWooslev fc Weaveilll995h 
from massive stars and so build up quickly during ac- 

^ The data presented herein were obtained at the W.M. Keck 
Observatory, which is operated as a scientific partnership among 
the California Institute of Technology, the University of California 
and the National Aeronautics and Space Administration. The Ob- 
servatory was made possible by the generous financial support of 
the W.M. Keck Foundation. 



five star formation epochs, while Fe-peak elements are 
produced in all supernovae. [a/Fe] abundance patterns 
are therefore particularly useful and have been central to 
developing our current picture of the assembly and star 
formation history of the Milky Way and spiral galaxies 
in general fsee lPritzl et al.ll2005l and references therein). 

Bright, young stars record the current gas phase 
abundance in a galaxy; to probe the earliest forma- 
tion times, one must target older, lower mass, fainter 
stars. It is only recently that individual red giant branch 
(RGB) stars in our nearest neighbor galaxies in the Lo- 
cal Group have been within reach of the high resolu- 
tion spectroscopy needed for detailed chemical abun- 
dance analysis. These stars in Local Group dwarf 
galaxies show a much greater range of abundance ra- 
tios at all metallicities compared to stars in the Milky 
Wav halo (Venn et all [mil iShetrone et al.ll200lL 1200^ 
iGeisler et al.ll2005t iTofstov et al.l 120091 ) . suggesting that 
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they have had a much more comphcated star formation 
history than the halo's progenitor(s). Detailed abun- 
dances beyond the Milky Way and its nearest remaining 
neighbors are needed to establish the broader patterns of 
star formation in galaxies of all masses. 

Unfortunatel y, at ~780 kpc from the Milky Way 
(|Hollandlll99l . even M31 is distant enough that older 
RGB stars are too faint (V~23) to obtain the required 
high signal-to-noise ratio (SNR'^ 60) and high spectral 
resolution (A/AA ~ 20,000) to measure detailed abun- 
dances in individual stars. Fortunately, globular clus- 
ters (GCs) can also be targeted. Unhke single stars, 
high resolution spectra can be obtained of unresolved 
GCs out to ~ 4 Mpc with current telescopes. GCs 
are bright enough (—10 < My < —6 mag) and have 
low enough velocity dispersions (2 — 20 km s^^) that 
even weak lines (~15 mA) can be detected in spectra of 
their integrated light. Detailed abundances have never 
been obtained from unresolved GCs because techniques 
have not existed to analyze them. We have developed a 
new method for analyzing high resolution integrated light 
(IL) spectra of single age, chemically homogeneous stel- 
lar population s to obtain detailed e lement abundances 
as described inlBernstein fc Mc Willia m ( 2002, 2005) and 
iMcWilliam fc BernsteinI (|2008D . her'eafter ';MB08," 

Our method has been developed and tested on a "train- 
ing set" of Milky Way and LMC GCs with well deter- 
mined properties from studies of individual stars. IL 
spectra of the training set GCs were obtained by scanning 
a 32x32 arcsec^ region of the cluster cores in the Milky 
Way clusters, and a 12 x 12 arcsec^ region in the LMC. 
Note that slit scanning is only necessary for the resolved 
GCs we target in our training set, not for unresolved, 
extragalactic GCs. The training set GCs were chosen to 
cover the range of metallicity, horizontal branch morphol- 
ogy, mass, velocity dispersion, and age available in the 
Milky Way and LMC systems. Using this training set, 
we have compared abundances determined with our IL 
method and abundances from the literature determined 
for individual RGB stars. Based on the Milky Way train- 
ing set, we estimate empirically that our [Fe/H] abun- 
dances are accurate to within < 0.1 dex, and all other 
elemen t ratios ([X/Fe]) accurate to within < 0.1 dex (see 
IMBOSI and S. Cameron et al. 2009, in preparation). We 
also derive approximate ages (>10 Gyr) for our entire 
Milky Way training set. Using the larger age range of 
the LMC clusters, which includes clusters as young as 
~10 Myr, we have found that our method can clearly 
distinguish clusters over a large range in age. Our accu- 
racies at the youngest ages are described in a forthcoming 
paper (J. Colucci et. al 2009, in preparation). 

In addition to the fact that GCs are among the old- 
est stellar populations in galaxies, there is ample evi- 
dence that GCs are a good tool for tracing both the 
early formation of galaxies themselves and star forma- 
tion throughout their histories. The fact that the num- 
ber of GCs in a ga laxy scales w ith the total luminosity 
of the galaxy (e.g. lHarrig|fT991f ) suggests that the GCs 
trace total star formation. Additionally, young massive 
clusters are seen forming in regions with high star for- 
mation rates, such as mergers of gas rich galaxies (e.g . 
ISchweizer fc Seitzeij[l99l IWhitmore fc Schweizeijll995[ ). 
indicating that GCs form in major episodes of star for- 



mation, throughout the lifetimes of galaxies. While the 
old age of GCs alone suggest that they record the earli- 
est star formation episodes in galaxies, there is addition- 
ally strong evidence that both blue (metal-poor) and red 
(metal-rich) sub-populations of GCs in normal galaxies 
correlate with the host galaxy's lumin osity and overall 
metallicitvfe.g. [Brodie fc Stradeil [2006f) . This suggests 
that both sub-populations are a record of the formation 
history of the galaxy, with blue clusters possibly tracing 
the earliest star formation in dark matter halos and red 
clusters tracing the later formation after the gas is more 
enriched. Finally, GCs are also relatively easy spectro- 
scopic targets to analyze because, to first order, they are 
simple stellar populations (SSPs), and can be approxi- 
mated with a single age and single metallicity. For all 
of these reasons, detailed abundance analysis of GC sys- 
tems is a powerful tool for understanding the formation 
of galaxies beyond the Milky Way. 

M31 is the closest large galaxy to the Milky Way and 
an ideal target for galaxy formation studies using GCs. 
Like the Milky Way, M31 has a large system of GCs, most 
of which are old (> 10 Gyrs), and has a bimodal metallic- 
ity distribution with peaks at [Fe/H] ~ —1.4 and [Fe/H]^ 
-0.5 and a inean o f [Fe/H] ~ -1-2 (iBarmbv et al.l 120001 : 
iPerrett et 311120021 ). iBarmbv et all (|2007D find that M31 
and Milky Way GCs are structurally similar, with sim- 
ilar mass-to-light (M/L) ratios, and fall on a single GC 
fundamental plane. However, there are notable differ- 
ences that have been observed between the two galax- 
ies. To begin with, M31's GC system is more than a 
factor of 2 larger than the Milky Wav's ijGalleti et al.l 
20041) ■ has massive but diffuse GCs (|Huxor et al.l 120051 : 



Mackev et al.ll2006[ ). and metal-poor and compact GCs 



at ve r y large projected ga lactocentric radii ([Martin et al.l 
I2006t iMackev et al.l [20071 ). There is evidence for young 
and/ or intermediate age GCs in M 31 unlike in the Milky 
Way ([Beaslev et al.l 120041 : iPuzia eValJ 12001 . and a sig- 
nificant population of GCs of all metallicities kinemat i- 
cally associated with the thin disk ([Morrison et al.|[2"004D . 
There is also evidence for some chemical differences in 
M31 GCs; compared to the Milky Way, M31 GCs show 
enhanced CN molecular lines fe.g. iBurstein et aLlll984l ). 
and the first estimates of [a/Fe] ratios in M31 GCs have 
indicated it may on average be ^0.1—0.2 dex lower tha n 
in the Milky Way ([Puzia et al.ll2005l : lBeaslev et al.ll2005( ). 
More detailed studies are required to understand what 
these similiarities and differences imply for the forma- 
tion history of M31. 

Much of what is known about M31 comes from low 
resol ution spectroscopi c methods like the "Lick" system 
(e.g. iFaber et al.lll985f l. which ha ve been used to target 
spatially unresolved GCs (e.g. Huchr a et al.lll99"l[ ) to ob- 
tain the constraints to date on [Fe/H| and [a/Fe] ratios 
(e.g. iPuzia eFall 120051 [Beaslev et al.ll2005[ ). These line 
index systems were originally developed for studying the 
integrated light of galaxi es, which have velocity disper- 
sions of 100-300 km s- ^Faber fc .Tacksonl[l976l ). With 
the low velocity dispersions of GCs, individual spectral 
lines can be resolved with only slightly greater line blend- 
ing than one finds in individual RGB stars. The accuracy 
of the line index systems depends on calibrations that 
are sensitive to abundance ratios and overall metallic- 
ity. These limitations make detailed abundance analysis 
from individual lines, as in our method, an important 
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TABLE 1 
M31 Clusters 



Cluster 


RA 


Dec 


V 


E(B-V)=' 


Mvtot" 




[Fe/H]° 


[a/FcY 


HBRS 




(2000) 


(2000) 








(kpc) 








G108-B045 


00 41 43.26 


41 34 21.8 


15.83 


0.10 


-8.95 


4.87 


-0.94±0.27 


0.22 dz 0.19 


0.14 


G219-B358 


00 43 18.01 


39 49 13.5 


15.12 


0.06 


-9.53 


19.78 


-I.83dz0.22 


0.00 dz 0.27 


0.78 


G315-B381 


00 46 06.47 


41 21 00.2 


15.76 


0.176 


-9.24 


8.8 


-I.22dz0.43 






G322-B386 


00 46 26.94 


42 01 52.9 


15.64 


0.13 


-9.24 


14.02 


-I.21dz0.38 


0.25 dz 0.22 


0.41 


G351-B405 


00 49 39.81 


41 35 29.4 


15.20 


0.08 


-9.52 


18.2 


-1.80dz0.31 




0.71 



References. — a. I Rich et al] ||2005D, b. iFan et al .' (2008), c. Reddening corrected and using a distance 
modulus for M 31 of (m-M)= 24.47 i Hollandlll998l ). d". Pcrrctt et al, (2002) e. Low resolution spectroscopic 
metallicities of iHuchra et al.l (^1991), f. Low resolution [a/Fe] of IPuzia et al.l (:2005), g. Horizontal branch 
morphology ratios (HBR), defined as the "Mironov Index" = B/B + R, where B and R corres pond to the 
number of horizont al branch stars bluer or redder than V — I = 0.5 in the observed CMD of iRich et al.l 
(|2005l ) Other data: iBarmbv et aTl (|2000l ) 



TABLE 2 

Observation Log and Estimated SNR 



Cluster 


Exposure (s) 


SNR (pixel- 
4400 A 6050 A 


7550 A 


G108-B045 


16,200 


40 


60 


80 


G219-B358 


10,500 


40 


60 


70 


G315-B381 


14,300 


30 


50 


60 


G322-B386 


12,600 


30 


50 


60 


G351-B405 


10,800 


40 


60 


80 



next step. 

In this paper we present the analysis of 5 GGs in M31 
using high resolution spectroscopy of their integrated 
light. These represent the first set of clusters which we 
have observed as part of an ongoing project to study the 
GG population in M31 with the goal of constraining the 
stellar populations and formation history of M31. From 
these clusters, we derive detailed stellar abundances of 
old populations in M31 for the first time. We obtain re- 
sults for 14 elements: Mg, Al, Si, Ga, Sc, Ti, V, Gr, Mn, 
Fe, Go, Ni, Y, and Ba. We also present ages, velocity dis- 
persions, and reddening constraints for this sample. In 
§ [21 we describe our observations and data reduction. In 
§[3] we present analysis of the velocity dispersions and im- 
plications for our equivalent width abundance analysis. 
In § |4] we describe our abundance analysis method and 
present abundance results in § [SJ In § [5] we discuss the 
M31 GG chemical abundances in the contexts of galaxy 
formation, low resolution spectra line index abundances, 
and overall consistency with broadband photometric col- 
ors and existing Hubble Space Telescope (HST) GMDs. 

2. OBSERVATIONS AND REDUGTIONS 

In selecting GG targets in M31, we have focused ini- 
tially on clusters which are spectroscopically confirmed, 
and are estimated from low resolution indexes to have 
abundances in the range -1.8 < [Fe/H] < -0.5 dex. At 
the time of selection, this was the range in which we were 
most confident of the accuracy of our analysis based on 
our previous work wi th the Milky Way training set clus- 
ters described above ()MB08l and S. Gameron et al. 2009). 
The targets were further selected to be relatively well 
studied, reasonably isolated and well outside of M31's 
disk to reduce confusion and the chance of confusion 



Fig. 1. — A single frame taken with the Keck I guide camera is 
shown here to illustrate the relative size of a GC in M31 (half-light 
radius ~1 arcsec) relative to the 1.7 X 7 arcsec^ slit. The guider 
images a reflection of the sky off the polished slit plate, so that the 
slit itself is clearly visible where no image of the sky is reflected. 

from interloping sources. We also avoided the bright- 
est, most massive GGs while still targeting GGs bright 
enough to get sufficient SNR in a few hours of obser- 
vations, as they will have the highest velocity disper- 
sions and thus broader, less pronounced spectral lines,. 
With this in min d, we have selected th ese GGs from the 
Barmby catalog (|Barmbv et al.ll200Clf ). This initial set 
of GGs has magnitudes between 15 and 16. Magnitudes 
and spatial information are listed for all of the GGs in 
Table [1] along with low resolution abundance estimates 
from the literature and horizontal branch morphologies 
from HST imaging when available. 

For this work we obtained high resolution IL spectra 
of the M31 globular clu sters using the HIRES spectro- 
graph (jVogt et al.lll994[ ) on the Keck I telescope over the 
dates 2006 September 10-14. We used the D3 decker, 
which provides a slit size of 1" .7 x 7.0" and spectral res- 
olution of i? = 24, 000. The GGs in thi s sample have 
half light radii (r/i) between ^0.6" — 1.1" (|Barmbv et al.l 
i 2007i). We ca lculate from surface brightness profiles in 
iBarmbv et all (|2002l . |2007[ ) that -70-90% of the fight 
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Fig. 2. — Example M31 GC spectra in Lick Mgb region shown decreasing in [Fe/H] from our analysis (top to bottom). Note that this 
region is dominated by saturated lines and is shown for illustration as a region familiar in low resolution spectra analyses. 



fell in the 1" .7 x 7.0" slit for each cluster. An illustration 
of the observing setup is shown in Figure [TJ The wave- 
length coverage of the HIRES spectra is approximately 
3800-8300 A. Exposure times were between 3—4 hours 
for each GC and are listed in Table [2] along with SNR es- 
timates at three wavelengths, one in each HIRES CCD. 
Data was reduced with the MAKEE software package 
available from T. Barlow^. In analysis of more recent 
observing runs, we have compared data analyzed with 
MAKEE and with the HiRes Redux pipeline produced 
by J. X. Prochaska, which has many routines and strate- 
gies in common with the MIKE Redux pipeline. While 
HiRes Redux produces lower noise spectra overall, and 
often traces weak orders more accurately, the MAKEE 
results are accurate and sufficient for our analysis, partic- 
ularly since we explicitly avoid regions with sky emission 
or absorption features entirely. Example spectra for the 
five M31 GCs are shown in Figure[2l where it can be seen 
that many individual Fe I, Fe II, Mg I and Ti II lines are 
easily identified. 



3. VELOCITY DISPERSIONS 

One of the strengths of our ILS analysis is the amount 
of information and the number of constraints available 
in the high resolution spectra. In addition to checks re- 
lated to the abundance measurements themselves (see 
§ 14. 3p . we also have overall photometric colors, mag- 
nitudes, and internal kinematics (velocity dispersions). 
Mass estimates from velocity dispersions can also help to 
constrain the basic initial mass function, contributing to 
the overall consistency of our understanding of the stel- 
lar population. Measurements of the velocity dispersions 
also tell us the spectral line resolutions we can expect 

^ MAKEE was developed by T. A. Barlow specifically for re- 
duction of Keck HIRES data. It is freely available on the W eb at 
I http://spider.ipac. caltech.edu/staff/tab/makee/index. html I . 



for our abundance analysis. We compare these results to 
measurements in the literature when available. 

We measure velocity dispersions from the high resolu- 
tion IL sp ectra of the M31 GCs f ollowing the method de- 
scribed in iTonry fc David (|1979( ) , which is implemented 
in the IRAF^ task fxcor. The GC spectrum is cross cor- 
related on an order by order basis with a suitable tem- 
plate star. The full width at half-maximum of the cross 
correlation peaks (FWHMcp) measured by fxcor is then 
converted to a line-of-sight velocity dispersion ( aots ) us- 
ing an empirical relation between the two. This relation 
is established by cross correlating the original template 
star spectrum with artificially broadened versions of it- 
self that are made by convolving it with Gaussian profiles 
corresponding to aobs of 2—25 km s~^. RGB stars are 
suitable template stars for old stellar populations; in this 
analysis, we use a spectrum of HR 8831 (type G8 III) 
taken during the 2006 September run as the template. 

In Table [3] we report aobs and the associated l-tr er- 
rors measured for orders between 4000—6800 A. Only 
orders with high SNR and weak atmospheric absorp- 
tion are used in the cross correlation. We also avoid 
orders that include t he saturated Balmer lines. Recently, 
IStrader et al.l (|2009D have found a weak trend of aobs with 
wavelength in a fraction of the GCs they observed. We do 
not find any correlations between aobs and wavelength for 
G108, G315, or G322. Ho\ Yever a small c orrelation exists 
for G219 and G351. Like IStrader eTall (f2009l we find 
aobs decreasing from blue to red orders by ~1 km s~^. It 
is unclear if this is due to a color/metallicity mismatch 
of the GCs and the template star (these two GCs are the 
more meta l-poor ones in our sa mple) or other systematic 
effect (see IStrader et al.l I2009L for a larger discussion) . 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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TABLE 3 

Velocity Dispersion Measurements 



Cluster 




RMS 


^ord 


Error 


<^best 


RMS 




Error 






(km s"-"-) 








(km s"-*^) 








(km s^^) 


G108-B045 


10.24 


0.44 


26 


0.09 


10.14 


0.37 


9 


0.12 


9.82 ± 0.18 


G219-B358 


10.99 


1.00 


13 


0.28 


10.36 


0.84 


5 


0.38 


8.11 ± 0.36 


G315-B381 


9.87 


0.43 


28 


0.08 


9.65 


0.41 


10 


0.13 


< lO'' 


G322-B386 


11.42 


0.58 


19 


0.13 


11.35 


0.37 


6 


0.15 


11.49 ± 0.24 


G351-B405 


12.29 


1.06 


16 


0.27 


12.01 


0.41 


7 


0.15 


8.57 ± 0.45 



References. — (a) iDiorgovski et al.1 (|1997l ). fbl iPetersoiil (|1989l 'l 

Note. — All measured velocity dispersions are given without application of aperture 
corrections (see §[3]). aobs is the mean measured from all usable orders, abest is the mean 
measured for higher SNR orders between 4800—5800 A. 



Because of this issue, in Table |3] we have reported (ibest , 
which is the dispersion obtained for a subset of the high- 
est SNR orders between 4800-5800 A that do not show 
correlations with velocity dispersion and also give the 
smallest RMS errors. The number of orders used in each 
measurement are recorded in columns 4 and 7 of Table 

m 

Four of the GCs observed here have aobs mea- 
sur ements in the lite rature from previous observations 
bv IDiorgovski et al.l ([19971 . The spectra used by 
IDiorgovski et al.l ( 19971 ) were taken with a slightly nar- 
rower slit (1".15 X 7.0"), which will lead to a velocity 
dispersion roughly 0.5 km s^^ larger, which is com- 
parable to the measurement errors. We find that our 
measurements for Gl() 8 and G322 agree with those in 
IDiorgovski et al.l (|1997D within the quoted errors, but our 
observed velocity dispersions for G219 and G351 are 2 — 4 
km s^^ larger. It is po ssible that the diff e rence between 
our results and those of IDiorgovski et all (|1997D are due 
to a difference in wsini of the template star used in the 
cross c orrelation. The template used bv lDjorgovski et all 
(fT997t ) has wsim= 7 km s ^. We measure a mean FWHM 
of 8 km s^^ from line widths in our template star. Note 
that this is a combination of the intrinsic stellar line 
width (including wsini) and the instrumental resolution, 
which is FWHM = 7.7 km s~"'^based on an analysis of 
arc lines taken through the appropriate slit. This sug- 
gests that the rotational velocity of the star is quite small 
(less than or equal to roughly 1 km when added in 
quadrature) . The lower wsini of our template could there- 
fore cause the higher inferre d apbs found here . Both our 
measurements and those of IDiorgovski et all (|1997D re- 
quire an additional correction of approximately -1-14% to 
convert apbs to a projected central cr^ as described in 
IDiorgovski et al.l ([1997). This is the geometrical correc- 
tion appropriate to obtain a systemic line of sight velocity 
dispersion from a velocity dispersion measured within a 
radius of 2 — 3 x r^,. 

For completeness, we note that iBarmbv et al.l (|2007f ) 
have predicted aperture velocity dispersions for these 
same GCs by modeling surface brightness profiles. They 
have compared their predictions to observations of veloc- 
ity dispersions in the literature to derive an empirical cor- 
rection between the two. We note that our velocity dis- 
persion measurements for G219 an d G351 are mo r e con - 
sistent with the trend deri ved bv IBarmbv et al.| (|2007f) 
than the measurements by IDiorgovski et al.n 1997f ). A 
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Fig. 3. — My-u relation for Milky Way GCs (open circles), with 
data from Pryor & Meylan (1993) and Harris (1996). Training set 
Milky Way GCs are highlighted in gray squares and the five M31 
GCs are shown as filled triangles. Training set and M31 GCs are 
plotted with reddening correct e d V m agnitudes and velocity dis- 
persions from lPrvor &: MevlanI 111993 ). The line width parameter 
R calculated for a range of velocity dispersions is shown on the 
inset axis. 

velocity disp ersion of 12.1 ± 1.3 km s^^ has also been 
measured by IStrader et al.l (j2OO90 , which is more consis- 
tent with the value we measure. 

Finally, we report the first velocity dispersion measure- 
ment from high SNR, high spectral resolution, IL spectra 
for G315-B381. Our measurement of 9.65 ± 0.13 km s~^ 
is consistent with the upper limit of 10 km s~^ found 
bylPetersonI ()1989f ). 

Our velocity dispersions, when combined with GC Vh 
measurements in the literature, allow us to make order of 
magnitude estimates of the cluster M/L ratios, r/i for all 
of the GGs in t his sample, wit h the e xception of G315, 
can be found in IBarmbv et al.l (j2007D . A simple calcu- 
lation, following iSpitzed (I1987I ). assuming the virial the- 
orem and an isotropic velocity distribution for the GGs 
results in dynamical M/L~1.5— 2.5 for these four GGs. 
Here we have assumed M = 2.5v^r/G, where r, the half- 
mass radius, is 1.3xr^, and v, the three-dimensional ve- 
locity dispersion is \/3 x a. We can also compare the dy- 
namical M/L t o those predict e d wit h population synthe- 
sis models. In iPercival et al.l (|2009r ). the Teramo group 
provides M/L ratios based on their isochrones and an ap- 
propriate IMF from iKroupal (|200lD . As we use the Ter- 
amo isochrones for our abundance analysis (see § 14. 2p , 
we can use isochrones which are self-consistent with our 
later analysis and well-matched to each cluster. Using 
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our own abundance and age constraints, the most ap- 
propriate isochrones (ages and metallicities) for our clus- 
ters have a M/L~1.9— 2.9 based on the Teramo group's 
population synthesis work. It is typical to obtain sys- 
tematically different M/L from dynamical and popula- 
tion synthesis techniques; this difference is probably due 
to the inclusion of low mass stars in the population syn- 
thesis estimate, which are actually ejected in GCs due to 
dynamical evolution. We find a ratio between dynami- 
cal and population synthesis M/L ratios of 0.84 ± 0.27, 
which is consistent with the value of 0.73 ± 0.25 found 
by Barmby ct al. (2007) for a l arger sample of M31 GCs. 
[McLaughlin fc van der Marell (|2005| ) find a similar ratio 
of 0.82 ± 0.07 for a sample of Milky Way and old LMC 
GCs. These M/L values further identify these four M31 
GCs as consistent with the familiar Milky Way GC pop- 
ulation. 

The relationship described above between light and ve- 
locity dispersio n is shown in Figu r e[3l T he aobs for Milky 
Way GCs from iPrvor fc MevlanI (1993"), along with ab- 
solute magnitudes {My) listed by Harris (1996) are plot- 
ted as open circles. For comparison, we also plot the 
redd ening-corrected rnagnitu des and velocity dispersions 
from lPrvor fc MevlanI (|1993f ) for our training set GCs as 
gray squares and those measured for the M31 GC sample 
as black triangles. From Figure [3] we see that the Uobs 
measured here are consistent with what we would expect 
given the absolute magnitudes of this set of GCs. 

The velocity dispersion that we measure above quanti- 
fies another important point about the clusters, which is 
the limiting resolution that wc can obtain for individual 
spectral lines due to velocity broadening. The velocity 
dispersions we measure for this set of M31 GCs of 9—12 
km s^^ give a line width parameter (i? = A/FWHM) 
R =14,000—10,000, from which it is clear that we are 
fully sampling the line profiles with the R = 24, 000 spec- 
tral resolution provided by the slit configuration used 
here (the HIRES D3 decker). However, with dispersions 
this large we expect to have difficulty measuring equiv- 
alent widths for elements with weak lines (<30 mA), as 
we discuss in § [3 

4. ABUNDANCE ANALYSIS 

Our new method for obtaining detailed abundances 
of from their integrated light was developed and tested 
using a training set of Milky Way and LMC globular 
clusters. The basic method was described in detail in 
iBernstein fc McWilliaml (|2002f ). lBernstein fc McWilliamI 
(|200,'Tfl . lMB0^ and the full training set is presented in S. 
Cameron et al. (2009) and J. Colucci et al. (2009). The 
method is briefly summarized below. 

4.1. Equivalent Widths and Line Lists 

We measure absorption line equivalent widths (EWs) 
for individual lines in the IL spectra using the semi- 
automated program GETJOB (McWilham ct al. 1995b), 
with which we fit low order polynomials to contin- 
uum regions and single Gaussian profiles to individ- 
ual lines and double or triple Gaussians to line blends 
when necessary. Line lists and oscillator strength s were 
taken f rom McWiUiam & Rich (1994), Mc Willia m et all 
(199^), iMcWilliam (1998), IMB08I and IJohnson et all 
( 20061 ). We measure fewer lines than in standard in- 
dividual RGB star analyses, as lines in IL spectra are 



broader and weaker than in individual RGB stars due 
to the velocity dispersions of the GCs and the presence 
of continuum flux from warm stars. The lines and EWs 
included in our flnal analysis are listed in Table |4l As ex- 
pected, we find fewer clean lines for the GCs with higher 
velocity dispersions. 

4.2. CMDs and EW Synthesis 

In order to synthesize IL EWs to compare to our 
observed IL EWs, we next need to model the pop- 
ulation using theoretical single age, single metallic- 
ity isochrones. During analyses of our training set 
GCs, we performed extensive testing of a variety of 
isochrones from the Padova' ^ (iGi rardi e t al.l [2000f) and 
Teram o^ (iPietrinferni et al.l [2004 . 2006: ICordier et all 
I2OOI groups (see IMB08D . Because we require a 
large, self-consistent parameter space of both scaled- 
solar and a— enhanced isochrones, we have chosen to 
use the isochrones from the Teramo group for our 
abundance analyses. The isochrones available cover 
abundances from Z= 0.000 1—0.04 for both scaled so- 
lar and a— enhanced ratios. We choose to use the 
recommended canonical evolutionary tracks including 
an extended asymptotic giant branch (AGB) and 
a— enhanced low— temperat u re op acities calculated ac- 
cording to [Ferguson et al.l (|2005D . We also choose 
isochrones with mass— loss parameter of 77=0.2 because 
comparison with our training set GCs (particularly those 
of intermediate metallicity) show that they more accu- 
rately match the CMD of GCs; isochrones with 7/=0.4 
over-predict the fraction of extreme blue horizontal 
branch (HB) stars at intermediate metallicities. Given 
that we find blue HB stars are less critical to accurate 
abundance analysis than red HB stars (see ? 16. 4p . it is 
more important to our analysis that the isochrones ac- 
curately reproduce the red horizontal branch than that 
they populate the blue region of the horizontal branch 
when it may be present. 

We apply an IMF accordin g to the multip le— part 
power— law form described in iKroupal (120021 ) . which 
changes index at O.SMq and O.O8M0. Because we 
only observed the core ('^0.1-0.2xr/j) regions in our 
training set GCs, we removed stars less massive than 
'--^O.7M0 from the IM F to match the present day core 
mass functions ()MB08l ) that have experienced dynamical 
mass segregation and evaporation of low mass stars (e.g. 
Haumgardt fc Makino 2003). In this set of M31 GCs we 
have observed regions corresponding to ~1.5-3r^, well be- 
yond the core region where significant mass segregation 
is expected. Although we do expect present day GCs 
to be stripped of stars less massive than ~O.3M0 due 
to dynamical evaporation, the Teramo isochrones stop 
at O.5M0. We do not believe that neglecting stars in 
the O.3-O.5M0 range is a problem for our analysis, as 
stars less massive than O.5M0 contribute only ~l-2% 
to the total fiux of the population and <1% in absorp- 
tion features. By combining the model isochrones with 
cluster-specific IMFs, we can create synthetic CMDs for 
the range of possible ages and metallicities for which we 
have isochrones. Each synthetic CMD is divided into 

* Padova isochrones downloadable at |http: //pleiadi.pd. astro. it /| 

^ Teramo isochron es downloadable at 

|http: / / albione.oa-teramo.inaf.i"t7[ 
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TABLE 4 

Line Parameters and Integrated-Light Equivalent Widths for M31 GCs 



Species 


A 

(A) 


E.P. 

(eV) 


log gf 


EW(mA) 
G108 


EW(mA) 
G322 


EW(mA) 
G315 


EW(mA) 
G351 


EW(mA) 
G219 


Mg I 


4571.102 


0.000 


-5.691 




120.0 




86.7 


34.5 


Mg I 


4703.003 


4.346 


-0.666 




139.2 




99.8 




Mg I 


5528.418 


4.346 


-0.341 








147.3 


58.2 


Mg I 


5528.418 


4.346 


-0.341 








143.6 




All 


3944.016 


0.000 


-0.638 










134.9 


Al I 


3944.016 


0.000 


-0.638 










140.1 


Si I 


7405.790 


5.610 


-0.660 


58.4 










Si I 


7415.958 


5.610 


-0.730 


66.4 




50.1 






Si I 


7423.509 


5.620 


-0.580 


84.3 


66.5 


82.4 







Note. — Table 4 is presented in its entirety at the end of the text. Lines listed twice 
correspond to those measured in adjacent orders with overlapping wavelength coverage. 



~25 boxes of stars with similar properties, with every 
box containing ~4% of the total ^-band flux. The prop- 
erties of a flux-weighted "average" star in each box are 
used for the atmospheric parameters needed in synthe- 
sizing IL flux- weighted EWs. 

Flux- weighted synthesized EWs of lines a re calculated 
using our routine ILABUNDS (see iMBOSf ) which pro- 
duces an integrated light EW composed of the ~25 repre- 
sentative stars in ea ch GMD using spectral synthesis rou- 
tines from MOOG ()SnedeiJfl973D and model stellar at- 



G108-B045 



mospheres from Kurucz (e.g. ICastelli fc Kurucd [20041 )^. 
The synthesized EWs of each of the ~25 representative 
stars are averaged together, weighted by their respec- 
tive contribution to the total flux of the cluster. The 
assumed abundance in the line synthesis is adjusted iter- 
atively until the synthetic flux-weighted EW matches the 
observed IL EW to 1%. Initial abundance calculations 
are performed using scaled— solar Teramo isochrones 
and Kurucz ODFNEW stellar atmospheres, and then 
recalculated with a-enhanced Teramo isochrones and 
AODFNEW atmospheres when abundance results im- 
ply enriched a— element ratios are present. All five M31 
GGs analyzed here were determined to be a-enhanced, 
and thus the abundances we report use a-enhanced 
isochrones and AODFNEW atmospheres in all cases. 

All abundances were calculated under the assumption 
of local thermodynamic equilibrium (LTE). In the case 
of aluminum we also discuss the non-LTE correction sug- 
gested by Baumuellcr & Gchrcn (1997) in 15.31 

Our IL method as implemented h ere em ploys a fixed 
microturbulence law, as described in lMBOSl Since we do 
not adjust microturbulence values for individual stars in 
the synthetic CMDs we must be careful of line satura- 
tion. For this reason we only report abundances from 
lines with EW strengths less than ^--^150 mA for all el- 
ements except Fe II. Our analysis indicates that abun- 
dances from Fe II lines with EWs over ~100 mA start to 
deviate significantly from the linear portion of the curve 
of growth. To remain in the linear regime we avoid Fe 
II lines with EWs >100 mA wherever possible. In some 
GGs, like G351-B405, the only clean Fe II lines we mea- 
sure have EWs >100 mA. Using these lines in those cases 



2.0 Gyre 
3.0 Gyre 13Gi 
5.0 Gyre 15 G 



B.9'Gyrs ■ ^ 




^ The models are available from R. 
|http : / /kurucz . harvard, edu /grids . html | 



L. Kurucz's Website at 



-2.5 -2.0 -1.5 -1.0 -0.5 0.0 

[Fe/H] Isochrone 

Fig. 4. — Fe I (circles) and Fe II (crosses) abundance solutions. 
The input [Fe/H] value of the isochrone (plotted on the x-axis) 
equals the output [Fe/H] value of our solution where circles lie on 
the dashed line. 



will result in Fe II abundances that may be slightly high. 
In this work we will refer to lines we do not analyze 
because of large EWs as "saturated" (>150 mA in gen- 
eral, or >100 mA for Fe II). We note that, in principle, 
abundance upper limits could be obtained for elements 
for which all lines are "saturated." However, due to the 
flux-weighting of EWs from stars of different types, this 
analysis requires special care and is not investigated fur- 
ther here. 

We have also calculated hyperfine splitting (hfs) abun- 
dances for Al, Sc, V, Mn, Go, and Ba. For lines with EWs 
>20— 30 mA, desaturation by hfs can significantly reduce 
th e deriv e d abundances. We u se the hfs line lists given 
in iMBOSl iJohnson et al.l (|2006D . and references therein. 
Typical hfs abundances corrections here for Al, Sc, V, 
JVIn, Go and Ba were -0.1, -0.01, -0.1, -0.4, -0.2, and 
—0.15, respectively. 

4.3. Finding the Best-Fitting CMD 

In the course of our work with the IVIilky Way and LJVIG 
training set GGs, we have explored a variety of strate- 
gies for identifying the best-fitting GIVID ~ the GIVID 
which provides abundances that are most consistent with 
those obtained from the spectral analysis of individual 
RGB stars in those clusters. Our first efforts to iden- 
tify the best-fitting GIVID focused on obtaining a consis- 
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G315-B381 




-2.5 -2.0 -1.5 -1.0 -0.5 0.0 

[Fe/H] Isochrone 



Fig. 5. — Same as Figure|4]for G315. 



G322-B386 




-2.5 -2.0 -1.5 -1.0 -0.5 0.0 

[Fe/H] Isochrone 



Fig. 6.— Same as Figure|4]for G322. 



G351-B405 




-po l 

-2.5 -2.0 -1.5 -1.0 -0.5 0.0 

[Fe/H] Isochrone 

Fig. 7.— Same as Figure|4]for G351. 



tent [Fe/H] solution from the Fe I and Fe II lines. This 
strate gy was discussed with regard to the a nalys is of 47 
Tuc in [Bernstein fc McWiniamr(l200i 120051) . and|MBQ8|. 
However analysis of the full training set showed that the 
Fe I and Fe II lines typically give a self-consistent solu- 
tion at [Fe/H] values that are frequently more metal-rich 
than those obtained from analysis of individual stars. 
There are several possible explanations for a difference 
in [Fe/H] from Fe I and Fe II lines, although a detailed 
study of this problem is beyond the scope of this pa- 
per. We simply note here that a difference between Fe 
I and Fe II abundances in individual stars due to non- 



G219-B358 




-:^n r .. I .... I .... I .... I .... I .... ' 

-2.5 -2.0 -1.5 -1.0 -0.5 0.0 

[Fe/H] Isochrone 

Fig. 8.— Same as FigureHfor G219. 

LTE overionization was noted by (|Kraft fc Ivans! 120031) . 
and that inaccuracies in Fe II oscillator strengths can 
lead to large uncertainties in abundances (see recent dis- 
cussion in[Mclcnd cz fc Barbuv 2009). These results have 
led us to focus on a different strategy. Using the training 
set spectra, we have found that the best-fitting CMD can 
be consistently identified by taking advantage of the fact 
that the metallicity dependence of RGB morphology is 
reasonably well understood (Gallart et al . ,2005) . After 
extensive testing, we have found that we obtain consis- 
tent, accurate abundances by requiring that the abun- 
dance used in calculating the isochrones themselves be 
consistent with the abundance recovered by our analysis 
for the Fe I lines. This is consistent with the fact that 
we find the RGB to have the dominant influence on the 
strength of the Fe I spect ral line s, more so than on the 
Fe II lines, as discussed in lMBOSl . To clarify our analysis 
methods, we describe below the procedure we follow for 
each GC. 

For each CMD we calculate a mean [Fe/H] abundance 
from all available Fe I and Fe II lines. In this data set we 
measure 30-80 Fe I hues and 2-10 Fe II lines per GC. Fe 
abundance results for all CMDs for each GC are plotted 
in Figures H through [51 Circles and crosses correspond 
to Fe I and Fe II mean solutions. The horizontal axis 
shows the [Fe/H] of the a— enhanced Teramo isochrones. 
CMDs of the same age are connected by colored lines. 
Note that a— enrichment affects the [M/H] value of the 
isochrone; for clarity, the value plotted on the x-axis is 
the true [F e/H] value rather than t he overall metallicity 
[M/H] Csee lPietrinferni et al.l[200fih . 

As mentioned above, our criteria for selecting a best- 
fitting CMD is that the Fe abundance calculated from 
the Fe I lines is consistent with [Fe/H] used to produce 
the isochrone. This criteria is met where the solutions 
cross the dotted black lines in Figures|3]through[51 When 
the Fe I solution that crosses the dotted line for a given 
age (color) lies between two isochrones, we interpolate 
an appropriate isochrone acco rding to the prescription 
recommended in lPietrinferni et al.l (j2006D . We then have 
7 possible solutions where isochrones intercept the dotted 
black line in Figure H] — one abundance solution for each 
age. 

We next isolate the best-fitting CMD out of these 7 us- 
ing diagnostics commonly used in standard stellar abun- 
dance analyses. These diagnostics concern the stability 
of the [Fe/H] solutions, which should not depend on the 
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Fig. 9. — Diagnostic plots for G108-B045, for ages of 15, 10, 7, 5, and 1 Gyr (top to bottom). Oldest solutions have smallest Fe I 
standard deviation and smallest dependence on EP, wavelength, or or reduced equivalent width ( log(EW/wavelength)) for this cluster. Fe 
I and Fe II lines are marked by dark circles and light squares, respectively. Gray points mark lines rejected by a sigma clipping routine 
when calculating the mean abundances. The solid line shows the linear fit to the Fe I lines and dashed lines show the 1 a deviation of 
points around the fit. Dark and light diamonds mark the final average Fe I and Fe II abundances. 



parameters of the individual line (c^xcitation potentials, 
wavelengths, or reduced EWs^). In principle, as in in- 
dividual RGB stars, these diagnostics reflect the accura- 
cies in the physical properties of the atmospheres used 
in the synthesis. The abundance versus excitation po- 
tential (EP) diagnostic is sensitive to the temperature of 
stars, while the abundance versus reduced EW diagnos- 
tic is sensitive to the microturbulcnt velocities of stars. 
The abundance versus wavelength diagnostic is poten- 

Reduced EW = log(EW / wavelength) 



tially sensitive to the age of the CMD, because stars of 
different temperatures dominate the IL flux at different 
wavelengths. Unlike in RGB stars, in an IL spectrum, 
correlations with EP and wavelength can be caused by 
an inaccurate temperature distribution of stars in the 
CMD, which, for example, could be the result inaccurate 
modeling of HB morphology in the isochrones. Likewise, 
correlations with reduced EW can be the result of inac- 
curate proportions of stars of different gravities as well 
as a symptom of an inaccurate microturbulent velocity 
law. These effects are difficult to unravel without ad- 
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Fig. 10. — Diagnostics for G315-B381. The smallest Fe I standard deviation and smallest dependence on EP, wavelength, and observed 
equivalent width at ages of 10-15 Gyr. Symbols are the same as in Figure [9l 



ditional constraints on the CMD, and identifying these 
is not the primary goal of this work. We therefore take 
the existence of any correlations merely as an indication 
that a given isochrone is less representative of the true 
CMD then one with weaker correlations. We attempt 
to identify a best-fitting CMD solution by selecting an 
isochrone that minimizes these correlations. We use lin- 
ear least squares fits between [Fe/H] and these parame- 
ters to identify and quantify the strength of any existing 
correlations. Plots illustrating the behavior of the Fe 
abundances with EP, wavelength, and reduced EW are 
shown in Figures [5] through [T31 From these plots, we 
obtain 5 diagnostics: the slope of [Fc/H] with EP, wave- 



length and reduced EW, and the standard deviation of 
the [Fe/H] solution for Fe I lines and Fe II lines. 

For any one of these diagnostics, there is not a statis- 
tically significant difference between the quality of the 
solution from CMDs within a range of ±5 Gyrs. For 
example, the slope of the [Fe/H] vs. EP relationship in 
Figure [TUl for G315 looks essentially the same for CMDS 
between ages 5 and 15 Gyrs. However, while the differ- 
ence in these diagnostics may be small over a wide range 
in CMD age, we do find that they change monotoni- 
cally, and are strongly correlated with each other. This 
suggests that there is clearly a preferred age and [Fe/H] 
range of CMD for each GC. 
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Fig, 11. — Diagnostics for G322-B386, The smallest Fe I standard deviation and smallest dependence on EP, wavelength, and observed 
equivalent width at ages of 7-13 Gyr. Symbols are the same as in Figure |9] 



To see this more clearly, we plot these diagnostics for 
the M31 sample in Figures [HI through [T51 These plots 
show all five diagnostics as a function CMD age for the 
7 CMDs that satisfy the original selection criteria (see 
Figures [4] through [8]) . From Figures [14] through [181 it 
is clear that for all of the GCs in our sample, all five 
diagnostics simultaneously imply that better solutions 
are obtained for CMDs with ages >7 Gyr. For these GCs, 
as for old Milky Way GCs analyzed as part of our training 
set (S. Cameron et al. 2009), we find the acceptable 
CMD ages typically cover a range of 5 Gyrs (e.g. 10—15 
Gyrs or 7—13 Gyrs). This is not surprising because the 
CMDs themselves change very little over those ranges in 



age. We discuss these age constraints in detail in the 
next section. 

From this range of acceptable ages, we select one CMD 
to use for a final analysis run of all elements for which 
we measure lines. For old GCs, such as the present sam- 
ple, our abundance results are quite insensitive to which 
CMD in this age range is used. Again, this is not surpris- 
ing as the CMDs in this age range are very similar. This 
weak dependence is quantified in Figure [TOl in which the 
upper plot shows the small difference in abundance be- 
tween the oldest and youngest CMDs in the acceptable 
age range for each GC (discussed in §[5TT]). Nearly all el- 
ements change by <0.1 dex, and older CMD ages always 
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Fig. 12. — Diagnostics for G219-B358. The smallest Fe I standard deviation and smallest dependence on EP, wavelength, and observed 
equivalent width at ages of 7-13 Gyr. Symbols are the same as in Figure |9] 



give smaller abundances. The lower plot of Figure [19] 
demonstrates that the derived abundance ratios are even 
more robust. Since the change in the abundances of most 
elements tracks the change in Fe, the net difference in the 
abundance ratios is <0.05 dex in almost all cases. 



5. RESULTS: CHEMICAL ABUNDANCES 

We have measured abundances from the available clean 
lines of a— elements, Fe and Fe— peak elements, and neu- 
tron capture elements for each GC. Final abundances, 
the number of analyzed spectral lines, linc-to-line scat- 
ter, and the age of the best-fitting isochrone are reported 



in Tables [5] through [51 All abundance ra tios relative to 
Fe use the solar abundance distribution of lAsplund et al.l 
(|2005f ). with a solar loge(Fe)= 7.50. Abundance ratios of 
Sc II, Ti II, Y II, and Ba II are reported with respect 
to [Fe/H]ii. Figures [201 through [Mj show the M31 abun- 
dance ratios (green circles) compared to our Milky Way 
training set IL abundances (red squares). Error bars for 
IL abundances in Figures [^Dj through [Mj correspond to 
the statistical error of the deviation in abundances from 
the Njines available for each species, as reported in Ta- 
bles [S] through [5] Note that these errors are often larger 
for the Milky Way training set abundances, which is a 
result of the smaller luminosity sampling of these GCs 





Eicitolion Potenliol (eV) Wovelerglh (A) loq( Eqjiuolenl WlOtH/Wovelenqlh) 

Fig. 13. — Diagnostics for G351-B405. The smallest Fe I standard deviation and smallest dependence on EP, wavelength, and observed 
equivalent width at ages of 10-15 Gyr. Symbols are the same as in Figure [9l 



(5-30 % of the total flux) and, in some cases, lower SNR 
spectra. 

5.1. Iron and Ages 

As in abundance analyses of individual stars, the first 
element we analyze is Fe because the large number of 
available transitions provide a wide variety of very useful 
diagnostics and consistency checks. As outlined in § 14.31 
we first use the Fe lines to constrain our best CMDs. We 
find this sample of M31 GCs is old, with preferred ages 
>7 Gyrs, which is consistent with age estimates from 
previ ous photornetric and low resolution spectroscopic 
work (jRich et al.ll2005l : lHuchra et al.lll991[ ). We also find 



this sample of GCs spans a metallicity range of [Fe/H]~ 
-0.9 to -2.2, which is within the range of [Fe/H] of Milky 
Way GCs system. Below, we briefly discuss the best age 
and [Fe/H] for the individual GCs in this sample. 

We flnd the preferred range of CMD ages for G108 to 
be 10-15 Gyrs. We pick a best CMD age of 15 Gyrs for 
our final abundance determinations. For the best solu- 
tions we find negligible trends of Fe I abundance with 
EP and a slight correlation remaining with wavelength 
and reduced EW, as shown in Figure [9l We note that 
the remaining correlation between Fe abundance and re- 
duced EW suggests that the microturbulent velocity law 
we have applied is not perfect for every star in the syn- 
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Fig. 14. — Fe line abundance diagnostics for G108. All diag- 
nostics are normalized to their maximum values so that they can 
be shown on the same scale. The thick black line is Fe I a, thick 
dashed orange line is Fe II a, dotted blue line is the slope in [Fe/H] 
vs. A, thin dashed green line is the slope in [Fe/H] vs. reduced 
EW, and the thin red solid line is the slope in [Fe/H] vs. EP. 



G315, Fel a = 0.172 




Isochrone Age (Gyrs) 

Fig. 15.— Same as Figure [141 for G315. 

G322, Fel a=0.150 




Isochrone Age (Gyrs) 

Fig. 16.— Same as Figure [141 for G322. 



thetic CIVID. However, while the correlation may still be 
present in the best solution, the overall trend in corre- 
lations over all ClVIDs is still very clear, and it is easy 
to select the most appropriate CIVID. It is important to 
note that the scatter in microturbulence values around 
the relation we have adopted is at least ±0.2 km s^^ 
for studie s of both dwarf and RGB stars in the M ilky 
Way (e.g. lBensbv et al.ll2005HFulbright et al.ll2006f ). We 
have experimented with adjusting the microturbulence 
law within this range, which could potentially reduce the 
Fe abundance standard deviation. However, as long as 




Isochrone Age (Gyrs) 

Fig. 17.— Same as Figure [Til for G351. 



G219, Fel tj = 0.203 




Isochrone Age (Gyrs) 

Fig. 18.— Same as Figure [Til for G219. 



the dependence of Fe abundance on reduced EW is small, 
this will not significantly alter the mean Fe abundance. 
We have not been able to find a microturbulence law that 
improves the overall abundance solution (i.e. improves 
all five diagnostics discussed in ii l4.3p . A more detailed in- 
vestigation of this issue is beyond the goals of this study. 
To preserve the self-consistency of our solutions between 
GCs, we do not alter the original microturbulent veloci- 
ties of the best- fitting CIVID. A detailed discussion of the 
small (<0.1 dex) systematic error between the IL spec- 
tra abundance analysis and that for individual stars is 
included in S. Cameron et al. (2009). For the purposes 
of this work, we avoid systematic error issues by concen- 
trating on relative comparisons between IL abundances 
determined in the same way for 1VI31 and JVIilky Way GCs 
in our training set. To summarize our results, the final 
abundance solution for G108 is [Fe/H]= -0.94 ± 0.03, 
where the uncertainty is the standard error in the abun- 
dances from all Fe lines {cr / ^/ Nunes — !)■ 

The preferred CIVID age range for G315 is also 10-15 
Gyrs. Four of the five diagnostics show best solutions 
at the oldest ages, therefore we use the 15 Gyr solution 
as our best CIVID. The 15 Gyr solution for G315 has 
negligible trends of Fe I abundance with EP, wavelength, 
and reduced EW, which can be seen in Figure [101 The 
final abundance for G315 is [Fe/H]= -1.17 ± 0.02. 

As can be seen in Figures [TT] and [TBI for G322 we find 
the EP correlation to be slightly stronger at the oldest 
ages, resulting in a slightly younger preferred age range 
of 7-13 Gyrs. We use a best CIVID age of 13 Gyrs, and 
find a negligible correlation with wavelength and reduced 
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Fig. 19. — Robustness of the abundances over the r^S Gyr range 
of acceptable isochrone ages. Top plot shows the difference in the 
loge(X) of older and younger solutions. Bottom plot shows the 
corresponding difference in abundance ratios. Dotted lines around 
zero mark a change of ± 0.1 dex. 

EW, but a slight correlation between Fe abundance and 
EP for this solution. The final abundance is [Fe/H] = 
-1.14 ±0.03. 

G219 is the only GC in this sample which appears to 
be slightly younger than the others. FigurefTSlshows that 
four out of five Fe line diagnostics are best for ages of 7- 
13 Gyrs. We pick a best CMD age of 10 Gyrs, which is in 
the middle of this preferred age range. Figure [T^ shows 
that this best solution — the 10 Gyr CMD — still shows 
slight correlations of Fe abundance with EP and reduced 
EW. We find an [Fe/H]=-2.21 ± 0.03 for G219, which 
makes it one of the most metal-poor GCs in the Local 
Group confirmed by high resolution spectra to date. 

The preferred CMD age for G351 is 10-15 Gyrs. We 
note that G351 has a larger Fe I standard deviation than 
any other GC in the sample for any solution, which can 
be seen in Figure [T3l We choose a best CMD age of 15 
Gyrs, and note that this solution has small correlations 
of Fe abundance with EP and wavelength, as well as a 
fairly significant correlation with reduced EW. We find 
[Fe/H]= -1.33 ±0.04. 

5.2. Alpha Elements 

As described in § [TJ [a/Fe] abundance ratios are a 
valuable tool for studying the star formation history of 
a galaxy, a— elements are produced primarily in SNII 
that occur on timescales of 1-20 million years, which cor- 
responds to the lifetimes of massive stars. While Fe- 
peak elements are produced in both type la (SNIa) and 
SNII, the SNIa contribution dominates on timescales of 



TABLE 5 
G108-B045 Abundances 



Species logioe(X) a Error [X/Fe]^ Nn^ 
15 Gyrs 



Si I 


7.15 


0.17 


0.12 


+0.58 


3 


Cal 


5.64 


0.17 


0.07 


+0.26 


7 


Sc II 


2.54 


0.26 


0.18 


-0.00 


3 


Ti I 


4.10 


0.21 


0.10 


+0.14 


5 


Ti II 


4.86 


0.34 


0.17 


+0.47 


5 


VI 


2.90 


0.24 


0.17 


-0.16 


3 


Or I 


4.82 






+0.12 


1 


Mnl 


4.05 


0.15 


0.09 


-0.40 


4 


Fe I 


6.56 


0.22 


0.03 


-0.94 


49 


Fe II 


6.99 


0.01 


0.01 


-0.51 


2 


Co I 


4.30 


0.19 


0.13 


+0.32 


3 


Nil 


5.47 


0.27 


0.14 


+0.18 


5 


Y II 


1.45 






-0.25 


1 



1 For Fe this quantity is [Fe/H]. 



TABLE 6 
G315-B381 Abundances 



Species 


logioe(X) 
15 Gyrs 


a 


Error 


[X/Fc]i 


N, - 

^ " lines 


Si I 


7.01 


0.27 


0.27 


+0.67 


2 


Cal 


5.51 


0.18 


0.07 


+0.37 


9 


Sell 


2.29 


0.38 


0.22 


+0.21 


4 


Ti I 


4.07 


0.15 


0.11 


+0.34 


3 


Ti II 


4.53 


0.30 


0.15 


+0.60 


5 


VI 


2.72 






-0.11 


1 


Cr I 


4.63 


0.15 


0.15 


+0.16 


2 


Mn I 


3.75 


0.14 


0.08 


-0.47 


4 


Fe I 


6.33 


0.17 


0.02 


-1.17 


61 


Fe II 


6.53 


0.50 


0.19 


-0.97 


6 


Nil 


4.97 


0.19 


0.09 


-0.09 


6 


YII 


1.44 






+0.20 


1 


Ba II 


1.45 






+0.25 


1 



1 For Fe this quantity is [Fe/H]. 



TABLE 7 
G322-B386 Abundances 



Species 


logioe(X) 
13 Gyrs 


a 


Error 


[X/Fc]i 


N, - 

^ " lines 


Mgl 


6.70 


0.08 


0.08 


+0.35 


2 


Si I 


6.70 






+0.27^ 


1 


Cal 


5.46 


0.16 


0.06 


+0.33 


8 


Ti I 


4.21 


0.14 


0.10 


+0.49 


3 


Ti II 


4.33 


0.21 


0.12 


+0.51 


4 


VI 


3.02 






+0.20 


1 


Cr I 


4.68 






-0.04^ 


1 


Mnl 


3.61 


0.07 


0.05 


-0.60 


3 


Fe I 


6.36 


0.16 


0.03 


-1.14 


35 


Fe II 


6.43 


0.41 


0.21 


-1.07 


3 


Ni I 


5.15 


0.26 


0.19 


+0.10 


3 


Ba II 


1.63 






+0.54 


1 



1 For Fe this quantity is [Fe/H]. 

^ A correction of —0.1 dex has been applied (see 

§[52J 

^ A correction of —0.25 dex has been applied (see 
§[53} 
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TABLE 8 
G219-B358 Abundances 



Species logioe(X) cr Error [X/Fe]^ Nn. 
10 Gyrs 

Mgl 

Al I 

Ca I 
Sc II 
Ti II 

Cr I 

Fe I 
Fe II 
Ba II 



1 For Fe this 

^ A non-LTE correction of +0.6 dex has been ap- 
phed. 

>109 years (e.g. ISmecker-Hane fc Wvsill992[ ). Thus, at 
early times many a— elements are produced while total 
Fe abundances are low, resulting in an enhanced [a/Fe] 
abundance ratio at low [Fe/H]. After the onset of SNIa, 
the total Fe abundance increases at a faster r ate than 
that of a— elements, decreasing the [a/Fe] ratio (jTinslevI 
Il979f) . Most GCs in the Milky Way have [a/Fe] ratios 
that are enhanced with respect to solar abundance ratios, 
similar to Milky Way halo stars at comparable [Fe/H]. 
This implies that Milky Way GCs formed when the ISM 
was dominated by enrichment by SNII. Like Milky Way 
GCs, we find all the GCs in our M31 sample to be en- 
hanced in Ca, Ti and Si. 

Abundances for Ca I come from 7—9 clean lines per 
GC, with rms scatter about the mean of 0.1—0.2 dex, 
which is similar to the scatter in our Fe abundances. 
We measure 4-10 Ti I and Ti II lines per GC, with 
slightly higher line-to-line scatter that may be due to 
weak blends. We are able to confidently measure 1—3 Si 
I lines in three of the GCs; lines in G219 and G351 were 
too noisy or badly blended to use. The one Si I line we 
measure in G322 is partially blended. To estimate the 
effect of weak blends on the Si abundance from the EW, 
we have used the SYNTH routine from MOOG, which we 
have modified to synthesize IL spectra. We estimate that 
at most our EW abundance measurement is '^O.l dex too 
high after visual inspection of synthesized IL spectra of 
different abundances. 

In Milky Way field stars, the [Mg/Fe] ratios behave 
similarly to Ca, Ti and Si. However, in Milky Way GCs, 
[Mg/Fe] sh ows inter- and intra- cluster abundance vari- 
ations (see iGratton et al.l I2004D with respect to other 
a— elements and star-to-star differences within individ- 
ual GCs. Analysis of our full training set of Milky Way 
GCs has revealed lower [Mg/Fe] than [Ca/Fe], [Ti/Fe], or 
[Si/ Fe] in t hree out of the six GCs where it was measured 
(see lMBOSi and Cameron et al. 2009). 

We have been able to measure [Mg/Fe] in three of the 
GCs in M31. Similar to the Milky Way training set GCs, 
in M31 we measure [Mg/Fe] to be lower ( [Mg/Fe] < 
dex) than other a— elements within individual GCs in 
two out of the three GCs where it was measured. In the 
other two clusters the Mg I lines had strengths > 150 
mA, and were therefore not analyzed in this work. 

We have performed spectral synthesis tests of the Mg I 
lines to see if the [Mg/Fe] depletion can be explained by 



TABLE 9 
G351-B405 Abundances 



Species logioe(X) cr Error [X/Fc]^ Nu. 
15 Gyrs 



0.13 +0.01 4 

0.06 +0.32 9 

-0.34 1 

0.23 +0.37 4 

0.25 +0.62 2 

-0.13^ 1 

-0.57 1 

0.04 -1.33 42 

0.10 -1.04 2 

-0.2^ 1 

0.20 +0.26 3 



^ For Fe this quantity is [Fe/H]. 

^ A correction of —0.25 dex has been applied (see 

§[M| 

3 Estimate from spectral synthesis of 6767 A (see 
Figure ^ 



TABLE 10 

Mean IL a Abundances for Milky Way 
Training Set and Current Sample of M31 
GCs 





Milky Way 


M31 


[Ca/Fe] 


+0.35 ± 0.08 


+0.34 ± 0.05 


[Ti/Fe] 


+0.46 ± 0.15 


+0.47 ± 0.10 


[Si/Fe] 


+0.52 ± 0.20 


+0.51 ± 0.21 


[Mg/Fe] 


+0.18 ± 0.39 


+0.15 ± 0.17 


[oCaTiSiMg/Fe] 


+0.38 ± 0.15 


+0.37 ± 0.16 


[aCaTiSi/Fe] 


+0.44 ± 0.09 


+0.44 ± 0.09 



line-to-line measurement error. An example of this test 
is shown in Figure [251 for the unblended 5528 A Mg I 
hue in the most metal-poor GC G219. From the EW of 
this line, we measure an abundance of [Mg/Fe]= +0.01. 
Overplotted are synthesized spectra with abundances of 
[Mg/Fe]= +0.3, 0.0, -0.3. G219 is extremely metal poor, 
and no other elements contribute to the line EW. It is 
clear from Figurc[25]that the closest matching abundance 
is [Mg/Fe] = 0.0, and that the line is inconsistent with the 
[a/Fe]= +0.3 measured for Ca and Ti, clearly showing 
that measurement error cannot explain the Mg deviation 
from the other a-element abundances in this cluster. In 
§ 16.21 we further discuss light element variations in the 
IL of GCs, and in it 16. 3. 21 we discuss further implications 
for low resolution IL abundances. 

We can address the star formation history of M31 by 
calculating an average [a/Fe] ra tio for each GC simi- 
lar to that in iPritzl et al.l (|2005f ). While these authors 
use Ca, Ti, Si, and Mg in their average [a/Fe], we 
note that [Mg/Fe] is probably not a good [a/Fe] indi- 
cator in the IL of GCs, for the reasons discussed above. 
Therefore, for the mean [a/Fe] for each GC discussed 
here, we include only Ca I, Ti I, Ti II, and Si I. The 
mean [a/Fe]= +0.36 ± 0.20, +0.50 ± 0.16, +0.40 ± 0.12, 
+0.44 + 0.16, and +0.49 ±0.13 for G108, G315, G322, 
G351 and G219, respectively, which is significantly and 
consistently enhanced relative to solar in all five M31 
GCs. We can also calculate the mean ratio for the indi- 



5.42 


0.14 


0.14 


+0.09 


2 


4.24 


0.05 


0.05 


+0.69^ 


2 


4.49 


0.14 


0.05 


+0.39 


8 


0.95 






+0.05 


1 


3.33 


0.22 


0.13 


+0.58 


4 


3.28 


0.29 


0.12 


-0.15 


7 


5.29 


0.21 


0.03 


-2.21 


47 


5.35 


0.12 


0.05 


-2.15 


4 


0.04 


0.08 


0.04 


+0.02 


4 



lan+it-ir ia We /HI 



Mgl 


6.22 


0.22 


Cal 


5.30 


0.16 


Sc II 


1.67 




Ti I 


3.94 


0.40 


Ti II 


4.49 


0.25 


Cr I 


4.44 




Mnl 


3.49 




Pel 


6.17 


0.26 


Fe II 


6.46 


0.17 


Ni I 






Ba II 


1.39 


?Q 
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Fig. 20. — a— element ratios. Gray points show Milky Way Stars and black points show mean Milky Way GCs abundances from single 
stars. Data for single stars is from Pifenn et afl JSOOJ), [Pritzl ct al. (2005) and refer ences ther e in. W hen possible, these abundance ratios 
have been adjusted to be consistent with the solar abundance distribution of .Asplund et aLl (120051) that was used in our analysis. Red 
squares show the abundances from our Milky Way training set ILS analysis and green circles show M31 ILS abundances. The average of 
[Ti/Fe]i and [Ti/Fc]n are plotted vs. [Fe/H]i, where we have measured them. 



vidual a— elements across the sample of GCs, and com- 
pare these values to similar means in our sample of Milky 
Way training set GCs. This comparison of abundances 
derived only with our IL spectra method avoids any po- 
tential sources of systematic error. The mean values for 
Ca, Ti, Si, and Mg are presented in Table fTOl In addition, 
we present the mean and deviation of all the [a/Fe] ratios 
including and excluding [Mg/Fe]. Table [TOl shows that 
GCs in both the Milky Way and M31 have extremely con- 
sistent [a/Fe] ratios. GCs in both galaxies are also con- 
sistent with t he Milky Way ha lo average a-enhancement 
of -1-0.35 fe. g. iMcWiliiamll 1997f l . The obvious implication 
of the [a/Fe] abundances in this small sample of M31 
GCs is that M31 was (or its now- merged components 
were) dominated by enrichment by SNII when these GCs 
formed. 

5.3. Aluminum 

Al abundances have been used to put constraints on 
chemical evolution models of the Milky Way because Al 
enrichment is p articularly sensitive to the details of SNII 
explosions (see iGehren et al.l 120061 ). Al abundances for 
Milky Way disk stars with [Fe/H]< —2 have some contri- 
butio n from explosive burning in SNII (|Andrievskv et al] 
120081 ). so that hke [a/Fe], Al a bundances reach a plateau 
value of [Al/Fe] - +0.3 fe.g. lBensbv et al.ll2005| ). How- 
ever, some stars in Milky Way GCs are found to be even 
more enriched, reaching levels as high as [Al/Fe]'^ +1 



(jCratton et al.ll2001h . Also, hke Mg, Al exhibits inter- 
and intra-cluster variations, which we discuss further in 
§ I6.2( suggesting that the influences on Al abundance are 
more complicated in GCs than in the field. 

We have only been able to make one measurement of 
Al I in this first sample of M31 GCs. This measurement 
was made from the 3944 A Al I line in the metal-poor 
GC G219. We were able to make two independent mea- 
surements of this 3944 A line because it was present at 
the ends of two adjacent orders. Both measurements give 
a consistent result of [Al/Fe]= -|-0.09±0.05. However, we 
note that Al I abundances deri ved from the 3944 A reso- 
nance line can be problematic: iMc William et a l. (1995^) 
find that the line is significantly blended with CH lines in 
some giant stars, which ca uses derived abundance s to be 
too high, and according to ' Baumueller fc Gchrcn' (' 19971 ) 
Al I abundances from this line will b e un derest imated 
by ^0 .6 dex due to non-LTE effects. lAndrievskv et al.l 
(|2008l ) find that this correction may be even larger in 
metal-poor hot stars. We do not see evidence in the IL 
spectrum for contamination of the 3944 A line by signif- 
icant CH blends. An appropriate non-LTE correction 
of +0.6 dex raises the abundance to [Al/Fe]= +0.69. 
This is higher than the [Al/Fe] ~0 halo average at an 
[Fe/H]=— 2.2, but consistent with the significantly en- 
hanced Al I found in some individual GC stars. This 
enhanced [Al/Fe] is also consistent with Al abundances 
we derive for three Milky Way GCs in our training set. 
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Fig. 21. — Aluminum abundances. Symbols are the same as 
in Fig ure 1201 GC data is from references compiled in ICarrettal 
ll200gi . A non-LTE correction of 0.6 dex has been added to G219 
abundance. 

The training set [Al/Fe] measurements provide much 
stronger evidence that the enhanced [Al/Fe] we mea- 
sure in GC IL spectra is real because they are not mea- 
sured from the potentially problematic 3944 A Al I fea- 
ture. The training set [Al/Fe] (see Figure [2T|) were 
measured fr'om the 6696/6698 A Al I doublet, which 
should not be contaminated by blends. In addition, 
iBaumueller fc GehrenI (|1997f) find that non-LTE effects 
are much smaller for these lines. The 6696 A feature 
is not detectable in G219, even though we measure an 
abundance as high as [Al/Fe]= -1-0.7 with the non-LTE 
correction. We performed spectral synthesis tests to 
check that the lower limit of [Al/Fe]= -1-0.7 is consistent 
with a 6696 A line that would be too weak to detect in 
the IL spectra. We indeed find that an [Al/Fe]= -1-0.7 is 
not high enough for the feature to be seen after convo- 
lution with the velocity broadening of G219 because it 
results in a line depth of < 0.01% of the continuum level. 

The 6696/6698 A features were also too weak in the 
other M31 GC spectra to get reliable abundance mea- 
surements from EWs, with the exception of G108, in 
which it unfortunately falls too near the end of an order 
for a good measurement to be made. The 3944 A line 
was saturated in all of the more mctal-rich GCs. 

5.4. Fe-peak Elements 

Ni, Cr, Mn Co, Sc, and V abundances are interesting 
because their p roduction generally tracks that of Fe, (e.g. 
Ilwamoto et al.l 1999), resulting in [X/Fe]^ 0. Therefore, 
deviations from [X/Fe]~ are particularly interesting 
because they imply special conditions may have been 
present, such as variations in the mass function of super- 
novae, variations of mixing of supernovae ejecta into the 
local ISM, metallicity dependent or explosion dependent 
supernova yields, or contributions from special types of 
supernovae (e.g. iMcWilliam 1997). 

Ni abundances in Milky Way field and GC stars gen- 
erally follow the expected Fe-peak element trend of 
[Ni/Fe]= at all metallicities. We are able to measure 
[Ni/Fe] in three M31 GCs, and find it to be consistent 
with the Milky Way abundance trend. Ni was measured 
from 3-6 Ni I lines in G108, G315, and G322. The three 
GCs have a mean <[Ni/Fe]>= -1-0.06 ± 0.14, which is es- 
sentially identical to the mean of the Milky Way training 
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Fig. 22. — Abundances fo r Fe peak elements Cr, Ni and Sc. Sym - 
bols arc the same as Figure [20] GC data from 'Kraft ct al.' (199^), 
[gnedcn ct al l il99%. Cohen ct al. (1999), Ivans ct al. (199!^), 
I^hetrone & Keand ('2000), Ivans ct al. (2001), Ramirez & Cohen 
(2002), Lec & Carney (2002), Sncdcn et al. (2004), Carrctta ct ah 
(12004), iCohenI (|2004 ). jTautvaisiene et al.l IQOO£), and iLee et al. 

mm 

set GCs, which have <[Ni/Fe]>= -H0.05± 0.21. 

In the most metal-poor cluster G219, all Ni I lines are 
too weak for EWs to be measured reliably. In G315 most 
Ni I lines are too weak or noisy for EW measurements, 
and the 7393 A line is blended with telluric absorption 
lines. We estimate from spectral synthesis tests of the 

noisy 6767 A line in G351, that [Ni/Fe] 0.2. The 

spectral synthesis and spectrum are shown in Figure 
where the uncertainty in abundance due to the noisy con- 
tinuum can be fully appreciated. 
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[Cr/Fe]= in Milky Way stars for [Fe/H]> -2, 
but [Cr/Fe]< for [Fc/H]< -2. The deviation from 
[Cr/Fe]= at fow metallicity may be due to different 
chemical enrichment for the lowest m etallicity stars in 
the Milky Way halo ()McWiman]|fl997l ). We measure Cr 
I abundances for the four most metal-rich M31 GCs that 
are consistent with the solar [Cr/Fe] average in both our 
Milky Way training set GCs and in Milky Way stars 
at these metallicities. Few [Cr/Fe] measurements exist 
for GCs with [Fe/H]< -2. When observed, [Cr/Fe] in 
GCs also follows the decreasing [Cr/Fe] a bundance trend 
observed in Milky Way halo stars (see iShetrone et al.l 
(HqOI) for M92 and NGC 2419, and lLetarte et al.l (|2OO60 
for M15 and Fornax GCs). We are able to measure 



Cr I in the low-metallicity M31 cluster G219, and find 
[Cr/Fe]= —0.15, which is also consistent with the de- 
creasing halo abundance trend below [Fe/H]^ —2. Cr 
abundances were calculated from one Cr I line in G108, 
G322 and G351, two Cr I lines in G315, and seven Cr I 
lines in G219. Most other Cr I lines in the more metal- 
rich GCs have EWs over 150 mA and were not analyzed. 
The Cr feature at 5409 A, which is the only Cr line we 
measure in G322 and G351, is partially blended with 
weak Ti I and Fe I lines, so that the EW abundance may 
be slightly high. We used spectral synthesis tests of the 
Ti I and Fe I blends around the 5409 A Cr I feature to 
estimate the effect of the blends on the derived Cr abun- 
dance. We find that the Cr I abundance derived with the 
original EW measurement may be approximately ^0.25 
dex too high. A correction of —0.25 dex to our EW abun- 
dances of [Cr/Fe]= -1-0.21 and -1-0.12, would result in 
[Cr/Fe]= -0.04, and -0.13 for G322 and G351, respec- 
tively. These [Cr/Fe] are consistent with the [Cr/Fe] 
of Milky Way GCs. 

Mn is a particularly interesting Fe— peak element be- 
cause unlike most Fe— peak elements the [Mn/Fe] ratio 
is not solar at most metallicities. From [Fe/H]= —1 to 
-2.5 dex, [Mn/Fe]-- -0.4, and increases to [Mn/Fe]= 
at [Fe/H]= 0. This trend is similar but opposite to that 
of a— elements (e.g. McWilliam 1997). Possible explana- 
tions for this trend are: Mn is underproduced in SNII 
and overproduced in SNIa, or that SNIa Mn production 
is dependent o n the metallicity of the progenitor star 
(|Grattonlll989( ). Recent observations of [Mn/Fc] in stars 
in the Sagittarius Dwarf Galaxy have shown that only the 
latter explanation can simultaneously r eproduce the data 
in bo t h the Milky Way and Sagittarius (|McWilliam et all 
I2OO3I: iCescutti et alll2008h . This demonstrates the im- 
portance of obtaining abundance ratios in a variety of 
environments fo r a f ull understanding of chemical evo- 
lution. Because iBergemannI ()2008h find that Mn abun- 
dances may be underestimated due to non-LTE effects 
over most of this metallicity range, we focus on a rela- 
tive comparison of our abundances with others calculated 
under similar assumption of LTE. 

We find that the Mn I abundances in the M31 GCs 
are consistent with our training set abundances and 
the Milky Way [Mn/Fe] abundance trend. Our abun- 
dances for Mn are measured from 3—4 Mn I lines in each 
the three most metal-rich GCs, G108, G315 and G322. 
[Mn/Fe] is measured from the Mn I 6021 A line only 
in G351. We do not see evidence for blends that would 
change the result by more than ~0.1 dex. 

The mean [Sc/Fe] in Milky Way stars and GCs is ap- 
proximately solar. We measure a mean [Sc/Fe] for 
GCs in M31, with a larger scatter between lines for in- 
dividual GCs and between the five GCs {a ~ 0.2 dex) 
than we see for other, easier-to-measure Fe— peak ele- 
ments (cr ~ 0.1). We find a similar mean and scatter 
in our training set GC abundances. We measure abun- 
dances for 1—4 Sc II lines in each of the five GCs. 

Co abundances in stars and GCs in the Milky Way 
track that of Fe for [Fe/H]> —2, so that over this range 
in metallicity the [Co/Fe]~ 0. We measure Co I from 
three lines for the most metal-rich GC G108. Co I fea- 
tures in the other four GCs are too weak to measure 
reliably, with the exception of the feature at 4121 A, 
which we find to be significantly blended. The G108 
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[Co/Fe]= +0.32 ± 0.13 is in agreement with Milky Way 
training set abundances at this metaUicity. 

Milky Way stellar and GC V abundances typically 
track the abundance of Fe, resulting in [V/Fe]~ 0. We 
measure a mean V abundance from G108, G315, and 
G322 of [V/Fe]= -0.02, which is consistent with our 
measurements of [V/Fe] in the Milky Way training set 
GCs. Abundances for G108 come from 3 V I lines, and 
abundances for G315 and G322 each come from the V I 
6081 A line. These V I lines are all weak (EWs-30-40 
mA), but it appears unlikely that there are any blends 
that would cause misleading abundance measurements. 

In summary, the Fe-peak element abundances in this 
sample of M31 GCs are similar to those measured in 
Milky Way GCs, and consistent with enrichment domi- 
nated by SNII. 

5.5. Neutron Capture Elements 

The relative abundances of neutron capture elements 
are important because they are particularly sensitive to 
the details of star formation, but also because the nu- 
cleosynthetic sources and yields of the rapid (r— process) 
and slow (s— p rocess) neutron capture reactions re main 
uncertain (see I Venn et all [200l iMcWiliimll [l99l . In 
particular, the difference between the Ba and Y abun- 
dances for stars in the Milky Way and stars in nearby 
dwarf spheroidal galaxies provides strong evidence for 
differences between the star formation histories o f large 
galaxies and their satellites (jShetrone et al.l [20011 120031 : 
iGeisler et all 120051 ). Perhaps more interesting is that 
these abundance differences are at times inconsistent 
with simple nucleosynthetic explanations, and thus can 
provide new cons traints on uncertain reaction sites (e.g. 
iVenn et al.ll200l . 

We are able to measure abundances for the strong lines 
of Ba II and Y II in some of the M31 GCs. Unfortunately, 
we are unable to measure abundances for weaker Eu II 
and La II lines due to the high velocity dispersions of 
this sample of GCs (see § [S]) . Typical Ba abundances in 
Milky Way GC stars are between [Ba/Fe]~ - 0.5 for 
[Fe/H]> -1 and [Ba/Fe]< for [Fe/H]< -1. Our Ba 
abundances for the M31 GCs are consistent with these 
trends and with what we find for our training set GCs; we 
also measure the lowest [Ba/Fe] for the lowest metallicity 
cluster G219. Ba II abundances come from 1-4 lines in 
each of the GCs except for the most metal-rich cluster 
G108, for which all line strengths were over 150 mA. 

For [Fe/H]> -1.5, Milky Way GCs typically have 
[Y/Fe]~ 0. We measure a mean value of [Y/Fe]~ for 
G108 and G315 in M31, as well, which is also consistent 
with what we find for the Milky Way training set. These 
Y abundances come from the 4883 A Y II feature. Be- 
cause the Y abundances are derived from a single line, we 
have performed spectral synthesis tests to confirm that 
the 4883 A line is relatively unaffected by blends. 

6. DISCUSSION 

Although this is only the first set in a larger sample 
of GCs from our M31 study, we already have several in- 
teresting results. First, in 16.11 we discuss the chemical 
enrichment history of the present sample of M31 GCs 
and compare it to the Milky Way and dwarf galaxy GC 
systems. In § 16.21 and § 16.31 we discuss the implications 



1.0 




-1 nl ...... T. .. I .............. i 

-2.5 -2.0 -1.5 -1.0 -0.5 0.0 
[Fe/H] 

1.01 I 



0.5 - 




-1 nl ......... i .............. i .. . 

-2.5 -2.0 -1.5 -1.0 -0.5 0.0 
[Fe/H] 

Fig. 24. — Abundances of neutron c aptu re elements. Symbols 
and references are the same as in Figure [20l 





1.10 








1.05 


. G219-B358 






1.00 






X 

3 






ized i 


0.95 






"o 

P 


0.90 






nori 


0.85 


: [Mg/Fe] = + 0.3 






0.80 


L [Mg/Fe]= 0.0 

: [Mg/Fe] = -0.3 


Mg I ■ 




0.75 







5524 5526 5528 5530 5532 

Fig. 25. — Spectrum synthesis test for Mg I 5528 A line in G219. 
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of our measurements for both GC formation and evolu- 
tion and IL abundance work at low resolution. In the 
final two sections wc comment on constraints that can 
be put on horizontal branch morphology and reddening 
of unresolved GCs using high resolution IL spectra. 

6.1. Chemical History of M31 CCs 
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TABLE 11 

Metallicity Comparisons 





ILS 


1 


2 


3 


4 


Cluster 


[Fe/H] 


[Fe/H] 


[Fe/H] 


[Fe/H] 


[Fc/H] 


G108-B045 


-0.94 ± 0.03 


-0.85 


-0.94 ± 0.27 


-1.05 ± 0.25 


-0.71 ± 0.11 


G219-B358 


-2.21 ± 0.03 


-1.92 


-1.83 ± 0.22 




-2.00 ± 0.11 


G315-B381 


-1.17 ± 0.02 




-1.22 ± 0.43 






G322-B386 


-1.14 ± 0.02 


-1.09 


-1.21 ± 0.38 


-1.62 ± 0.14 


-1.03 ± 0.22 


G351-B405 


-1.33 ± 0.04 


-1.77 


-1.80 ± 0.31 
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Fig. 26.— Spectrum synthesis test for Ni I 6767 A line in G351. 
Smoothed data is shown in red, and ovorplotted synthesized spectra 
from top to bottom correspond to [Ni/Fe]= -0.4,-0.2,0.0. The 
spectrum is noisy and continuum placement uncertain, but [Ni/Fe] 
is approximately ~ —0.2 with an uncertainty of ~0.2 dex. 



Overall, these five 1VI31 GCs are old and have chemical 
properties similar to those of most IVIilky Way GCs. All 
five GCs are enhanced in the a— elements Ca, Ti, and Si 
to the same extent as IVIilky Way GCs. Fe— peak element 
ratios are consistent with IVIilky Way abundance trends. 
These results are consistent with exi sting simple G alactic 
chemical enrichment scenarios fe.g. lTinslevlll979l) , in so 
far as it suggests that the gas in both the IVIilky Way and 
1VI31 halos was dominated by enrichment by SNII when 
these GCs formed. The similar levels of a— enhancement 
in this small sample so far suggests that 1VI31 and the 
Ivlilky Way are likely to have had similar IlVIFs and star 
formation rates at early times. Since we do not see ev- 
idence for the low [a/Fe] observed in GCs in the LIVIC 
or the disrupting Sagittarius Dwarf Galaxy, it does not 
seem likely that any of these five GCs are associated with 
recent satellite accretion events. 

6.2. Variation of Light Elements 

Variations in the abundances of light elements between 
individual stars within GCs have been observed in all 
GCs studied in detail since the phenomenon w as dis- 
covere d in the IVIilky Way GCs M3 and 1VI13 by ICohenI 
(11978). Of the light elements, Mg and O are observed to 
be depleted in some GC stars, while Na and Al are over- 
abundant. These abundance variations are related and 



have since been called the Na— O and Ivlg— Al anticorre- 
lations (see Gratton ct al. 20(34). Abundances varying in 
this way are predicted from high temperature (T > 10^) 
C— N— O cycle H— burni i ig ([Denisenkov fc Denisenkoval 
119901: iDenissenkov et al.l \19m . While these reaction 
products can in principle be brought up to the stellar 
surfaces during "deep mixing" on the RGB, the obser- 
vation of the abundance variations even in GC main se- 
quence stars has suggested that they are instead the re- 
sult of pollution by GC in termediate— mass AGB stars , 
for example, as discussed in lVentura &: D'Antonal (|2008l ). 
although this solution is not universal l y acc epted (see 
modeling of NGC 6752 in lFenner et"an (|2004f )). 

Star-to-star Mg variations have been more difficult to 
meas ure than those of O, Al, and Na in many IVI ilky Way 
GCs (jCarretta et al.ll2004D . iSneden et all (|2004D success- 
fully measured [Mg/Fe] variations in a large sample of 
stars in the IVIilky Way GCs 1VI3 and 1VI13, with values 
that range from []VIg/Fe]= —0.2 to -1-0.4. This variation 
is detectable in the IL as the mean []VIg/Fe] will be lower 
than the mean [a/Fe], as shown in § 15.21 We can also 
expect that our IL measurements correspond to a kind 
of average []VIg/Fe] for each GC, which is complicated by 
the extent of the Mg depletion present within the GC, 
and the fiux weighting of stars of different types. As a 
simple consistency check, we note that our IL [Mg/Fe] 
measurements fall within the range []V[g/ Fe]= —0.2 to 
+0A that is expected for individual stars (jSneden et al.l 
(2(304i )). The scatter in IL [Ivlg/Fe] between GCs tends 
to be high in both the IVIilky Way and 1VI31, and does not 
appear to correlate with any property of the GCs. These 
factors suggest that it is difficult to predict an expected 
value of []VIg/Fe] in the IL spectra of a GC at this time. 

AGB star pollution in GCs can also affect the abun- 
dance of [Al/Fe]. Abundances in some GC stars are ob- 
served to be as high [Al/Fe] > -1-1, significantly higher 
than the abundances seen in halo field stars. We de- 
termine an abundance for G219 of [Al/Fc]= +0.69 (see 
§ 15. 3p . which is close to the high value of [Al/Fe] ~ 1 
measured in some GC stars. Three of our training set 
IVIilky Way GCs have [Al/Fe] > -1-0.6, which was mea- 
sured from the more reliable Al 6696 A line, so we be- 
lieve that we are seeing evidence of the ]VIg— Al anticor- 
relation in IL measurements of [Al/Fe] in both the IVIilky 
Way and 1V131. 

Unfortunately, we were unable to measure oxygen 
abundances in this sample of GCs. This was either due to 
blends with telluric absorption lines or because the lines 
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Fig. 27. — Compariso n of high reso l ution ILS [Fo/H] with line 
ind ex measurements by IPerrett et al.l II2002I ) ( P02, o pen circles), 
[Huchra et al. (199j) (HB91, open squares) and lPuzia et al. (200|) 
(PPB05, open triangles) plotted as a function of V magnitude. 



were too weak to reliably measure from EWs with the 
reasonably high velocity dispersions of these GCs. Na 
lines were either too weak or very saturated. Future IL 
light abundances for GCs with smaller velocity disper- 
sions should be more useful for investigating the Na— O 
anticorrelation. 

6.3. Comparisons to Lick Indexes 

Much progress has been made using low resolution Lick 
index systems to measure global properties of GC sys- 
tems (see review by Brodie & Strader 2006). Low resolu- 
tion metallicities have helped establish the general trends 
in GC populations in other galaxies, and their impor- 
tance for tracing galaxy formation and evolution. In this 
section, we compare previous measurements of Lick index 
metallicities for our sample of GCs in M31 with the goal 
of understanding where and why abundances from low 
and high resolution spectra may differ. The Lick system 
or similar methods are critical to studying extragalactic 
GCs because low resolution spectra analyses will always 
need to be applied to systems that are so distant that 
high resolution analyses are impossible. 

In making these comparisons we caution that the defi- 
nition of [Fe/H] is slightly ambiguo us. Line index met al - 
licities are typically calibrated to the lZinn fc WestJ (Il98l l 
[Fe/H] scale as established for Milky Way GCs^. This 
metallicity scale is advantageous because it covers the 
metallicity range spanned by Milky Way GCs and can be 
easily applied to distant objects, providing a consistent 
IL metallicity scale for both Milky Way and extragalac- 
tic GCs. However, by definition this scale is based on 
blends of lines from multiple elements. This limits the 
information that can be reliably determined for individ- 
ual element abundances, including Fe. Another difficulty 
in these comparisons is that the calibration is based on 
the Milky Way GC system, which has a very consistent 
abundance pattern and HBR — [Fe/H] relation that is 
not necessarily the same in GCs in other galaxies. Since 
the calibration is based on spectral regions with blends of 



* The lZinn &: WestI 119841) scale is based on the integrated photo- 
metric parameter Qsg and measurements of absorption from Ca H 
and K, CN, Fe, and the Mgb region from low resolution integrated 
spectra and was calibrated to early high resolution abundance re- 
sults. 



several elements, it may be less accurate if targets don't 
have Milky Way-like abundance ratios. 

6.3.1. "[Fe/H]" 

Because of the ambiguities described above, the com- 
parison of Lick index metallicities to the high resolu- 
tion [Fe/H] determined in our analysis is interesting both 
when the abundances agree and when they disagree. Our 
IL [Fe/H] results for all five GCs are summarized in Table 
11, along with metallicity estimates in the literature. 

A comparison of the high resolution [Fe/H] derived 
from our analysis and low resolution Lick index [Fe/H] 
estimates is shown in Figure [27l plotted as a function 
of V magnitude for convenience. High resolution [Fe/H] 
are plott ed as solid symbols a n d the line inde x measure- 
ments oflHuchra et all (|1991h . IPerrett et all (2002 ) and 
iPuzia et al.l (|2005D correspond to open squares, circles, 
and triangles, respectively. Figure [57] shows that the 
Lick index [Fe/H] estimates agree within the errors for 
the higher metallicity clusters G108, G315, and G322. 
However, bigger differences appear at the lowest metal- 
licities; our measurement of [Fe/H]= —2.21 ± 0.03 for 
G219 is ~ 0.2 — 0.4 dex lower than previous estimates, 
and our measurement of [Fe/H]= —1.33 ± 0.04 for G351 
is ~0.5 dex higher than low resolution results (see Ta- 
ble [IT]). 

Larger discrepancies at the lowest metallicities are ex- 
pected for several reasons. First, line index strengths 
change little for [Fe/H]< —1.6, so calibrations to the 
Lick sy stem at low metallicity are uncertain (e.g. see 
iPuzia e t al. 2002). In the case of G219, a difference be- 
tween [Fe/H]=-1.8 and [Fe/H]=-2.2 would be partic- 
ularly difficult to detect at low resolution, emphasizing 
the importance of high resolution measurements at these 
lowest metallicity GCs. The discrepancy between the 
high resolution and Lick index abundance for G351 is a 
little more difficult to unde rstand, esp e cially since the 
CMD metallicity estimate of iRich et al.l ()2005D is similar 
to the Lick index result. This discrepancy stands out in 
our sample, however the difference is still within ^1.3cr 
of the quoted error of Huchra et al. (1991), and our own 
rms error for this cluster is slightly larger than for others 
in our sample. As a simple reality check, a visual compar- 
ison of the G351 spectrum and our other 4 GC spectra 
(see Figure [5]) suggest that G219 is substantially more 
metal poor than G351, while line index measurements 
put both G351 and G219 at [Fe/H] 1.8. 

It is possible that the Lick index analysis is compli- 
cated by a combination of factors that include both low- 
metallicity degeneracies and poor modeling of blue hori- 
zontal branches (see discussion of the effect of HB mor- 
phology on our results in 16. 4p . We note G351 ha s been 
obser ved to have a bimodal horizontal branch (Rich et al.l 
|2005[ ). We also note that G351 is one of the GCs we find 
to have depleted [Mg/Fe], suggesting that abundance ra- 
tio calibration degeneracies may be a particular problem 
for this cluster in the sample. 

6.3.2. [a/Fe] 

Recent progress has been made in developing SSP IL 
spectra models with variable element rati os for compar- 
ison \ yith Lick index a bsorption feat ures ([Thomas et al.l 
120031: iLee fc Worthevl[2005,: Schiavoni 120071) . In particu- 
lar, the models of mas et al. r pool " have made the 
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Fig. 28. — Comparis on of [g/Fe] f rom h igh resolution ILS with 
hne index estimates by IPuzia et ahl l|2005l ) (PPB05, open circles) 
plotted as a function of V magnitude. Red crosses correspond to 
mean [a/Fe] measured from Ca I, Ti I, Ti II and Si I only. Crosses 
are plotted with a +0.05 offset in V for visibility. Blue triangles 
correspond to mean [a/Fe] from Ca I, Ti I, Ti II, Si I and Mg I. 

[a/Fe] estimates of a portion of M3 1 GCs possible, in - 
cluding three of those studied here (iPuzia et al.ll2005h . 
In this case, [a/Fe] is determined from a comparison 
with models of Lick indexes with Fe— dominated and 
Mg— dominated absorption feature s. A significant result 
of the study of (|Puzia et al.ll2005f ) is that GCs in M31 
with ages >8 Gyrs have a mean [a/Fe]= +0.18 ± 0.05 
with a dispersio n of 0.37 dex, whi ch is 0.1 ~ 0.2 dex 
lower than what 'Puzia et al.' (2005") find for Milky Way 
GCs. However, jjeaslcy et al. (2005) find that the dis- 
crepancy between [a/Fe] in M31 and Milky Way GCs 
may be an SSP model— dependent result. 

Since we have found that [Mg/Fe] is likely to be de- 
pleted compared to other a— elements within GCs due 
to AGB star self— pollution, we expect that [a/Fe] ratios 
in GCs determined from indexes with Mg— dominated 
absorption features could be lower than the [a/Fe] abun- 
dances that we determine from Ca, Si, or Ti. 

To test this, we use the mean [a/Fe] from Ca, Si and Ti 
lines for each GC ( discussed in ^ 15. 2| a nd compare to the 
[a/Fe] estimated bv lPuzia et al.l (I2005D from L i ck ind exes 
in Figure [28l Measurements bv IPuzia et al.l (|2005l ) are 
plotted as open circles, and our mean [a/Fe] excluding 
Mg are plotted as red crosses. We measure a system- 
atically higher value for [a/Fe] than Puzia ct al. (200^) 
obtain for [a/Fe] in all three GCs. The largest discrep- 
ancy is_Joi_G219, fo^which we find [acaSiTi/Fe]= -1-0.47 
and IPuzia eFal] (|2005D estimate [a/Fe]= 0.0. We find 
the discrepancy for G219 is reduced, but not resolved, 
if we include our [Mg/Fe] = -1-0.04 measurement in the 
mean [a/Fe], which is shown by the blue triangle in Fig- 
ure Hi 

We note that, in addition to the degeneracy in the Lick 
indexes at low abundances, the effect of different [a/Fe] 
ratios in SSP modeling at low abundances is also very 
weak, further obscuring the resolution of the Lick index 
me asurements at low abu nda nces, as already dis cussed 
bv iMaraston et all ()2003D and IPuzia eFaH (|2005f ). It is 
likely that this is the cause of the remaining discrepancy 
in the [a/Fe] estimate from low resolution for G219. 

While our present sample size of high resolution IL 
abundances in extragalactic GCs is still small, it already 
removes some of the discrepancies between M31 and 



Milky Way GCs from Lick indexes that have been dis- 
cussed in the literature (jPuzia et al.ll2005t feeaslev et al.l 
l2005HBrodie fc Stradeill2006D . We find that for this sam- 
ple of GCs the true [a/Fe] from Ca, Si, and Ti appear 
similar for the two galaxies, and that an accurate esti- 
mate of this value in GCs for interpretation of the chem- 
ical enrichment history of a galaxy must come from el- 
ements other than Mg due to the peculiarities of Mg 
abundances in GCs. Our results indicate that the unex- 
pected, low [a/Fe] ratios in metal-poor GCs may be an 
artifact of the uncertainties in line index systems at low 
metallicities. 

6.4. Horizontal Branch Morphology 

In general, the position of a star on the HB is a func- 
tion of metallicity; it is expected that metal-poor GCs 
will have bluer HB morphologies than GCs with higher 
abundances. The fact that a number of Milky Way GCs 
at the same metallicity are observed to have very differ- 
ent HB morphologies has led to the conclusion that at 
least one important "second— parameter" plays a role in 
HB morphology (see review bv lGallart et al.ll2005( l. SSP 
models that reproduce the HBR— [Fe/H] relationship of 
the Milky Way are difficult to establish because HB mor- 
phology is sensitive to a variety of factors (e.g. mass loss , 
age, helium abundance, a nd others) (e.g. lLee et al.lll995 
iRecio-Blanco et aLll2006() . Universal SSP models will be 
even more difficult to establish if the HBR— [Fe/H] rela- 
tionship is different in other galaxies. There is evidence 
that this might indeed be the case; observations of GCs 
in both M31 and Fornax suggest the HBR— [Fe/H] rela- 
tionsh ip is offset to lower metallicities tha n in Milky Way 
GCs f.Rich et al.ll2005l : iSmith et aLlflOOSh . 

Low resolution studies have found it important to 
consider the effect of HB morphology on IL spectra 
of unresolved GCs because the presence of old, hot 
stars on the blue HB can mimic light from young 
main sequence stars, resulting in young or intermediate 
age determinations for GCs that are actually old (e.g . 
'de Frcitas Pacheco fc Barbuvl fl995l : iBeaslev et al.l 120021 : 
Schiavon ct al. 20041). In the next section we discuss the 
effect of HB morphology on abundance and age determi- 
nations with our high resolution IL method. 

6.4.1. Effect on Age and [Fe/H] Determinations 

We have used our training set Milky Way GCs to per- 
form tests to assess the effect of inaccurate proportions 
of red and blue HB stars on our results. The Teramo 
group has produced isochrones with two different values 
for mass loss during stellar post-main sequence evolu- 
tion, a parameter which influences HB morphology. Of 
the two possible values, we have used isochrones with the 
less extreme mass loss parameter. This mass loss param- 
eter produces very blue [B — V < 0.1) HB stars for our 13 
and 15 Gyr CMDs at [Fe/H] < -1.8, but does not pro- 
duce bimodal HBs, as is common in many intermediate 
metallicity GCs. Because we are using the less extreme 
mass loss parameter, the fear is that the HB may not be 
as blue as appropriate in some cases. 

To test the effect of this potential error, we have added 
blue HB stars by hand into our synthetic CMDs to as- 
certain the consequences of underprcdicting the number 
of blue HB stars with our choice of isochrones. In these 
tests we conserve the total number of stars and total 
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Fig. 29. — M31 GC spectra near H(5 region, shown in order of increasing HBR from top to bottom. 



flux of the HB. We find tfiat even though blue HB stars 
can contribute 10 — 15% of the total flux of the popula- 
tion at wavelengths below 5000 A, most of the Fe I lines 
are found over the range 4500-7000 A, which is not as 
strongly influenced by the blue HB stars. Empirically, 
the addition of these stars into the CMD changes the 
derived Fe I abunda nces by < -1-0.05 dex. This is not 
unexpected, because IMBOSI showed that the Fe I EW 
strengths are dominated 80%) by the luminous, cooler 
stars on the RGB, AGB and red HB, with little effect 
from hot turn off or hot blue straggler stars, particularly 
at redder wavelengths. 

As an example of an extreme test case, we replaced all 
of the red HB stars with the same number of blue HB 
stars in the 10 Gyr best-fltt ing CMD for G21 9. We know 
from the HST CMDs of iRich et all (|2005l) that G219 
has both red HB stars and a large number of blue HB 
stars that are not entirely represented in the best-fitting 
CMD. Using this extreme CMD, we find that the Fe I 
abundance changes by <0.05 dex, the scatter in the Fe 
abundances increases, and the correlations of Fe abun- 
dance with EP, wavelength, and reduced EW get larger. 
From this extreme example, it is clear that not only do 
blue HB stars have little impact on the Fe I solution, 
but that accurate modeling of the red HB results in a 
much more self-consistent solution from all of the Fe I 
lines. Therefore, we can expect that having roughly the 
correct number of red HB stars in the synthetic CMDs 
will be more important for our abundance determina- 
tions than the correct number of blue HB stars. Given 
the small effect of blue HB stars on Fe I abundances, and 
given that our analysis only constrains the CMD age for 
old GCs to a 5 Gyr age range, we find that the accuracy 
of the blue HB morphology has a no significant (or even 
detectable) impact on our results. We conclude that the 
ages and metallicities that we derive for old GCs with 
our abundance analysis method will not be significantly 



affected by synthetic CMDs that have inaccurate blue 
HB morphologies. 

As a further consistency check, we can look for quali- 
tative information about the temperature distribution of 
the preferred CMD solutions in the Fe abundance vs. EP 
plots discussed in S 14.31 A symptom of too few hot, blue 
HB stars in the isochrone compared to the real GC would 
be increasing Fe abundances for lines with larger EP. For 
G108, G315, and G322 we find no significant trend in Fe 
abundance with EP, which implies that the distribution 
of stellar temperatures in the CMD solutions for these 
GCs are very accurate. 

For G219 we find decreasing Fe abundances with in- 
creasing EP, which suggests in this case that we actu- 
ally do have too many hot stars in the isochrone for this 
metal-poor GC. However, we also find higher Fe abun- 
dances at larger reduced EWs, which implies that we 
may have too many dwarf (low microturbulent velocity) 
stars in the isochrone being used in the analysis. Since 
the dwarf stars have higher temperatures than the gi- 
ant stars, for G219 the excess of hot dwarf stars may 
cause the observed trends. Note that the dependence of 
Fe abundance on reduced EW may be associated with 
a small error in the microturbulence for some stars, as 
discussed in § 15.11 However, the mean Fe solution is not 
strongly affected by these weak trends, and so we do not 
pursue this issue further. 

The only GC in the sample that shows the symptom of 
not enough hot stars in the Fe abundance vs. EP plot is 
G351, although it also shows the largest positive depen- 
dence of Fe abundance on reduced EW. The diagnostics 
suggest that there are both too many dwarf stars and 
not enough hot stars in the best CMD for G351. Also, 
the standard deviation for [Fe/H] is larger than for the 
other four GCs. While the solution for this cluster shows 
somewhat less self-consistency than for the other 4 clus- 
ters from the diagnostics discussed above, we emphasize 
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TABLE 12 
E{B — V) Comparisons 





E{B 


-V) 




Cluster 


ILS Limit 


1 


2 


G108-B045 


0.13 


0.18 


0.10 


G219-B358 


0.06 


0.06 


0.06 


G315-B381 


0.08 


0.17 




G322-B386 


0.15 


0.21 


0.13 


G351-B405 


0.06 


0.14 


0.08 



Ref ere nces. — (1 . ) iFan et aLl 
(|200g ) (2.) iRich et al.l $200^ ) 



that the overall scatter in the solution is still quite small, 
suggesting that our solution has not been dramatically 
affected by unavoidable problems (e.g. interloping stars 
or internal age or abundance variations) and that the sta- 
tistical uncertainty is a meaningful estimate of the overall 
accuracy of the analysis. Moreover, any difficulties are 
not likely to be due to the HB morphology. 

6.4.2. Consistency with Photometry 

We are fortunate that the details of GC HB morphol- 
ogy do not have a large effect on the abundances or 
ages derived from high resolution IL spectra using this 
method. However, we can see the effect of HB stars on 
the temperature sensitive Balmer lines. We are still test- 
ing a method to constrain the HB morphology of unre- 
solved GCs using the Balmer lines in high resolution IL 
spectra (J. Colucci et al. 2009, in preparation). For the 
purposes of this work, we simply check for qualitative 
consistency of the Balm er line profi l es in the IL spectra 
with the HST C MDs of iRich et all (|2005l). HBR ratios 
from IRich et all (l2005l) for G108, G322, G351, and G219 
are listed in Table [U where a value of zero would cor- 
respond to a GC with a strongly red HB, and a value 
of unity would correspond to a GC with a strongly blue 
HB. From the IL spectra in the region of HJ, shown in 
Figure [^Sl it is clear that G219 and G351 have more 
prominent Balmer line wings, and thus a larger contri- 
bution of flux from hot stars than G108, consistent with 
e xpectations . 

IRich et al.l (|2005l ) measure an HBR-[Fe/H] relation 
from HST CMDs for 18 M31 GCs that is offset to lower 
mctallicities than the HBR— [Fe/H] relation in Milky 
Way GCs. One explanation proposed is that the most 
metal-poor M31 GCs are ~ 1 — 2 Gyr younger than sim- 
ilar GCs in the Milky Way. We note that our lower pre- 
ferred CMD age range for G219 is consistent with this, 
although we are unable to put strong constraints on ab- 
solute age. 

6.5. Reddening Constraints from Broadband B — V 

Colors 

In principle, an interesting additional constraint on our 
potential CMD solutions can be derived by comparing to- 
tal broadband {B — V)o colors with existing photometry. 
This comparison could help us eliminate potential CMD 
solutions where the integrated colors of the CMDs are 
inconsistent with the photometry (see MB08). In prac- 
tice, however, we find that reddening estimates for GCs 
in M31 are too uncertain to provide constraints on CMDs 
of the GCs discussed here. Rather, our abundance anal- 
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Fig. 30. — Integrated {B — V)o colors calculated from synthetic 
CMDs are shown as a function of [Fe/H] for each age. Black points 
show the [Fe /H] and age solutions determined from this analysis for 
G108. The horizontal solid line corresponds to the observed B — V 
color for the GC from G alletj et al, , 12004 ). with an arrow and dot- 
ted line to show the E{B — V) correction due to Galactic reddening 
of Schlcgel et al. (1998J. (B — V)o calcul ated with r e dden ing and 
metallicity of iFan et all li2008h (F08) and IRich erall ||2005(1 (R05) 
are plotted as open squares and triangles, respectively. 
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Fig. 31.— Same as Figure [30l for G322. 
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Fig. 32.— Same as Figure [30l for G315. 



26 



COLUCCI ET AL. 




-2.5 -2.0 -1.5 -1.0 -0.5 

[Fe/H] 

Fig. 33.— Same as Figure l30l for G351. 
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Fig. 34.— Same as Figure l30l for G219. 

ysis constrains the CMD with enough fidehty that we 
can actuaUy put some hmits on the reddening for these 
individual M31 GCs. We note that GCs have previously 
been used in this wa y to probe reddening in galaxies (e.g. 
iBarmbv et ani2000f ). 

The colors for our synthetic CMDs using Teramo 
isochrones are shown in Figures [301 through [Mj The 
broadband {B — V)q colors calculated from our synthetic 
CMDs are plotted against the [Fe/H] adopted in the Ter- 
amo isochrones for each GC. The trend of {B — V)o with 
[Fe/H] for constant age is emphasized by the solid colored 
lines. Black points correspond to the CMD and [Fe/H] 
solutions for each age described in ii l4.3l for each GC. The 
black solid line in each figure corresponds to the observed 
B—V color for the G C from the Revised Bologna Catalog 
(jGalleti et all [Mil . Black arrows show the minimum 
reddening for the M31 hne-of-sight of E( B -V) = 0.06 
which is due to Galactic extinction alone (jSchlegel et al.l 

For each individu al GC we show the recent redden- 
ing determination bv lFan et al.l ()2008[ ). who have derived 
these from correlations between optical and infrared col- 
ors and metallicity using a large set of M31 GC low res- 
olution spectroscopic abun dances. The (B — V)o calcu- 
lated with the reddenings of lFan et al.l (|2008f ) are plotted 



as open squares, using the low resolution spectroscopic 
[Fe/H] solutions for the x-axis location. For G108, G322, 
G351 and G219 we also show the reddenin g adopted by 
the deep HST study of iRich et all ()2005f ). which is an 
avera ge of se veral previous measurements. Reddenings 
of Ric h et all ([2005|) are plotted as open triangles at the 
[Fe/H] determined from their HST CMDs. 

For each GC we can place an upper limit on the redden- 
ing that is consistent with the preferred CMDs for each 
GC discussed in § 15.11 This limit is calculated using the 
{B—V)q color of the youngest CMD in the preferr ed age 
range and the i?- ^ color from lGalleti et all (|2004l). The 
limits are li sted in Tab le [12] with the values of iFan et al.l 
(|2008i) and iRich et al.l |2005) for comparison. We note 
that synthetic CMDs created with the Teramo isochrones 
that employ a higher mass loss rate than those used here, 
and thus have bluer HBs overall, would result in slightly 
lower predicted {B — V)o. For the ages relevant here, the 
effect HB morphology is at most A{B — V)o ~ —0.03, 
and an indication of the uncertainty in our reddening 
limit. 

In general, we find that our preferred CMDs and the 
reddenings th ey imply agree be st with the [Fe/H] and 
reddenings of iRich et al.l ()2005[ ) , rather than with the 
reddenings inferred from low resolution spectroscopic 
metallici ties. F or G108, G322 and G351 the E{B~V) de- 
riv ed bvlFan et al. (200^) are higher than those adopted 
by iRich et al.l (|2005i ). but it is clear from Figures [30] 
through[34]that they are still consistent with a n old (>10 
Gyr) population at the metallicities applied bv lFan et al.l 
(g008). 

Although G219 has been verified to be ~10 Gyr or 
older by HST photometry, by low resolution spectral line 
indexes, and by high resolution spectra abundances in 
this work, it has B—V colors substantially too blue 
for this age. This inconsistency has been noted in the 
Bologna catalog. In Table [T2l we list the Galactic E{B — 
V) — 0.06 as the only reliable constraint. We also note 
that the difference between the observed and predicted 
B ~ V tor G219 is too large to be explained by a deficit 
of blue HB stars in the synthetic CMD alone, but may 
be reduced with some contribution from relati vely rare, 
but lu minous, post-AGB or UV-bright stars fsee lMoehleH 
[2OOII) . 

For G351, the B - V of the 10 Gyr CMD is con- 
sistent with E{B — V) = 0.0, lower than the Galactic 
E{B — V) = 0.06, which we again list as the most reli- 
able constraint in Table [T^ This difference in color may 
be at least partially explained by missing blue HB stars 
in the sy nthetic CMD as c ompared to the observed HST 
CMD in IRich et all (|2005[ ). 

G315 is the only GC discussed here for which there are 
no HST constraints available. We find an E{B — V) < 
0.08 is consist ent with o u r ana lysis. The larger value 
determined by I Fan et aTl ([2008') is hard to understand 
because their derived reddening implies a. {B — V)o that 
is inconsistent with an old population, even though the 
high and low resolution [Fe/H] are very similar. It is 
possible that the dis crepancy can ari se from inaccurate 
photometry used bv iFan et al.l ()2008D to determine the 
E{B-V). 

7. SUMMARY 
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We have applied a new method for detailed abundance 
analysis of high SNR, high resolution IL spectra of unre- 
solved GCs to a sample of five GCs in M31. From over 
60-100 resolved spectral lines in each cluster we have de- 
rived abundances from EWs for Mg, Al, Si, Ca, Sc, Ti, 
V, Cr, Mn, Fe, Co, Ni, Y and Ba. We have used our 
abundance analysis to put independent constraints on 
the ages and reddening of these M31 GCs, and used the 
high resolution IL spectra to measure velocity dispersions 
to high precision. 

We find these 5 M31 GCs to be similar to the Milky 
Way GCs system in several respects. First, they arc >1G 
Gyrs old, and span a range in [Fe/H] of —2.2 to —0.9. 
Second, their Ca, Si, and Ti abundances are enhanced to 
similar levels as Milky Way GCs, and suggesting that the 
gas reservoirs from which they formed were dominated by 
products of SNII at the time of their formation. Finally, 
the Fe— peak and the neutron capture abundance ratios 
studied here also follow Milky Way abmidance trends. 

We have confirmed that light element abundance vari- 
ations between stars within GCs can effect abundances 
derived from high resolution IL spectra for Mg and possi- 
bly Al. We suggest that part of the large scatter in [a/Fe] 



measurements of cxtragalactic GCs using low resolution 
line indexes may be due to the effects of Mg— dominated 
absorption features on line indexes. 

We have demonstrated that a significant number of 
quantitative constraints on galaxy and GC formation 
and evolution can be made for unresolved GCs using this 
new abundance analysis technique. We have shown for 
the first time that abundance ratios fundamental to un- 
derstanding galaxy formation can be obtained for other 
nearby, massive galaxies. While we have intentionally 
targeted "typical" GCs for this first sample in M31, a 
larger selection of GCs is crucial for a complete picture 
of the GC system and formation history of M31. Future 
work is needed to investigate other possible differences 
in the M31 GC system compared to the Milky Way, i.e. 
young GCs, high metallicity disk GCs, and enhanced ni- 
trogen abundances. 
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M31 INTEGRATED LIGHT ABUNDANCES 



TABLE 13 

Line Parameters and Integrated Light Equivalent Widths for M31 

GCs 



Species 


A 

(A) 


E.P. 

(eV) 


log gf 


EW(mA) 
G108 


EW(mA) 
G322 


EW(mA) 
G315 


EW(mA) 
G351 


EW(mA) 
G219 


iVig i 




u.uuu 


-o.oyi 




izu.u 




00. ( 


"iA ^ 
04. 


Mgl 


4703.003 


4.346 


-0.666 




139.2 




99.8 




Mgl 


5528.418 


4.346 


-0.341 








147.3 


58.2 


Mgl 


5528.418 


4.346 


-0.341 








143.6 




Al I 


3944.016 


0.000 


-0.638 










134.9 


Al I 


3944.016 


0.000 


-0.638 










140.1 


Si I 


7405.790 


5.610 


-0.660 


58.4 










Si I 


7415.958 


5.610 


-0.730 


66.4 




50.1 






Si I 


7423.509 


5.620 


-0.580 


84.3 


66.5 


82.4 






Cal 


4318.659 


1.899 


-0.295 




109.6 








Ca I 


4425.444 


1.879 


-0.358 


127.9 








46.6 


Ca I 


5581.979 


2.523 


-0.555 


99.0 


77.6 


86.3 


64.2 




Ca I 


5588.764 


2.526 


0.358 




138.0 


137.6 


122.8 


59.8 


oa 1 


ooyu. izD 


9 t^9T 
Z , Z 1 


n c^V"! 

-U.O 1 1 




yD.D 


QO a 

oy .D 


TA 9 




Ca I 


5601.286 


2.526 


-0.690 


112.3 


91.9 


89.7 


53.3 


21.3 


Ca I 


5857.459 


2.933 


0.240 




124.0 


145.0 


111.9 


43.1 


Ca I 


6122.226 


1.886 


-0.320 








149.8 


69.0 


Ca I 


6162.180 


1.899 


-0.090 










89.0 


v_^a 1 
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1 1 A9 
-1. ±^Z 
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71 n 

1 i.U 






Ca I 


6439.083 


2.526 


0.390 








147.8 


57.9 


Ca I 


6471.662 


2.526 


-0.686 


98.9 




78.8 


70.0 




Cal 


6493.781 


2.521 


-0.109 






124.8 


94.3 


42.7 


Cal 


6572.795 


0.000 


-4.310 


82.5 










Cal 


7148.150 


2.709 


0.137 




133.5 


146.4 


131.3 




Sc II 


4246.837 


0.315 


0.240 






141.7 




78.3 


Sc II 


4670.413 


1.357 


-0.580 


84.3 


86.4 


86.9 
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Sc II 


6604.600 


1.357 


-1.480 


44.6 










Ti I 


4991.072 


0.836 


0.380 








134.2 




Ti I 


4999.510 


0.826 


0.250 


143.5 


122.0 


134.8 


129.7 




Ti I 


5039.964 


0.021 


-1.130 


115.9 










Ti I 


5210.392 


0.048 


-0.884 




140.1 








Ti I 


5401.379 


0.818 


-2.890 






9.4 






Ti I 


5648.565 


2.495 


-0.260 


17.3 










Ti I 


5866.451 


1.067 


-0.840 


88.4 


71.7 




46.3 
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- i .DoU 
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Ti II 
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Ti II 
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-2.460 


85.2 




73.6 






Ti II 


4501.278 


1.116 


-0.760 




150.1 








Ti II 


4563.766 


1.221 


-0.960 




155.6 


136.4 


126.3 


68.9 


Ti II 


4571.982 


1.572 


-0.530 








156.5 


76.7 


Ti II 


4589.953 


1.237 


-1.790 


82.6 


91.7 


77.2 




44.5 


Ti II 


5381.010 


1.566 


-2.080 


92.0 


60.6 


75.3 






V I 


6039.730 


1.060 
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39.5 
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4274.806 
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-0.231 
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n nnn 
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«7 1 


Cr I 


5206.044 


0.941 


0.019 










84.5 


Cr I 


5208.432 


0.941 


0.158 










87.2 


Cr I 


5208.432 


0.941 


0.158 










93.3 


Cr I 


5409.799 


1.030 


-0.720 




144.7 


134.0 


118.7 


36.8 


CrI 


7400.188 


2.900 


-0.111 


85.9 




76.9 






Mn I 


4754.039 


2.282 


-0.086 


103.9 


73.1 


80.4 






Mn I 


6013.520 


3.070 


-0.250 


63.7 


,38.1 


38.7 






Mn I 


6016.620 


3.070 


-0.216 


73.4 




44.7 






Mn I 


6021.820 


3.070 


0.034 


65.1 


44.5 


62.5 


36.6 




Fc I 


3878.027 


0.958 


-0.896 










85.3 


Fc I 


3899.719 


0.087 


-1.515 










143.7 


Fe I 


4063.605 


1.557 


0.062 
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TABLE 13 — Continued 



Species 


A 


E.P. 


log gf 


EW(mA) 


EW(mA) 


EW(mA) 


EW(mA) 


EW(mA) 




(A) 


(cV) 




G108 


G322 


G315 


G351 


G219 


Fe I 


4071.749 


1.608 


-0.008 










120.5 


Fe I 


4114.451 


2.831 


-1.303 


72.4 




50.7 






Fe I 


4132.067 


1.608 


-0.675 










93.4 


Fc I 


4132.908 


2.845 


-1.005 










45.3 


Fe I 


4147.675 


1.485 


-2.071 




109.8 








Fe I 


4154.505 


2.831 


-0.688 










33.6 


Fe I 


4156.806 


2.831 


-0.808 






131.7 






Fe I 


4157.788 


3.417 


-0.403 


82.9 




71.8 






Fe I 


4174.917 


0.915 


-2.938 




100.8 


75.8 




45.1 


Fe I 


4174.917 


0.915 


-2.938 






113.5 






Fe I 


4175.643 


2.845 


-0.827 


105.6 


102.9 


90.2 




20.0 


Fe I 


4181.764 


2.831 


-0.371 










50.3 


Fe I 


4182.387 


3.017 


-1.180 


82.1 




74.0 






Fe I 


4187.047 


2.449 


-0.514 




136.1 






56.1 


Fe I 


4191.437 


2.469 


-0.666 










59.5 


Fe I 


4195.340 


3.332 


-0.492 








123.5 




Fe I 


4199.105 


3.047 


0.156 






135.9 






Fc I 


4202.040 


1.485 


-0.689 










96.3 


Fc I 


4206.702 


0.052 


-3.960 


140.1 




115.6 


87.1 




Fc I 


4216.191 


0.000 


-3.357 








93.3 


47.6 


Fe I 


4222.221 


2.449 


-0.914 


120.7 


110.4 


119.6 




50.2 


Fe I 


4227.440 


3.332 


0.266 










73.7 


Fe I 


4233.612 


2.482 


-0.579 






142.9 






Fe I 


4250.130 


2.469 


-0.380 










69.4 


Fe I 


4250.797 


1.557 


-0.713 










86.4 


Fe I 


4260.486 


2.399 


0.077 










94.9 


Fe I 


4271.164 


2.449 


-0.337 










78.3 


Fe I 


4271.774 


1.485 


-0.173 










138.1 


Fe I 


4282.412 


2.176 


-0.779 








131.8 


86.8 


Fe I 


4325.775 


1.608 


0.006 










139.1 


Fe I 


4337.055 


1.557 


-1.704 




127.5 








Fe I 


4369.779 


3.047 


-0.803 




88.9 








Fe I 


4404.761 


1.557 


-0.147 










142.2 


Fe I 


4415.135 


1.608 


-0.621 










121.6 


Fe I 


4427.317 


0.052 


-2.924 










96.0 


Fc I 


4430.622 


2.223 


-1.728 






138.1 






Fe I 


4442.349 


2.198 


-1.228 




135.2 


144.2 






Fe I 


4443.201 


3.071 


-1.043 










21.7 


Fe I 


4447.728 


2.223 


-1.339 




114.4 


119.3 






Fe I 


4466.562 


0.110 


-0.600 




128.0 


139.1 






Fe I 


4494.573 


2.198 


-1.143 




141.3 








Fe I 


4602.949 


1.485 


-2.208 


142.7 


105.5 


113.7 


97.8 


45.4 


Fe I 


4632.918 


1.608 


-2.901 


107.6 




98.3 






Fc I 


4691.420 


2.990 


-1.523 




146.1 


76.1 


78.6 




Fe I 


4736.783 


3.211 


-0.752 


106.3 


107.2 


113.8 






Fe I 


4871.325 


2.865 


-0.362 




157.9 








Fel 


4872.144 


2.882 


-0.567 




158.2 




128.7 




Fe I 


4890.763 


2.875 


-0.394 








154.0 




Fel 


4891.502 


2.851 


-0.111 








141.9 




Fel 


4903.316 


2.882 


-0.926 




130.6 


137.8 


123.8 




Fe I 


4918.998 


2.865 


-0.342 








141.7 


60.9 


Fe I 


4920.514 


2.832 


0.068 










90.6 


Fe I 


4966.095 


3.332 


-0.871 


125.4 




113.8 


84.8 


78.8 


Fe I 


4994.138 


0.915 


-2.969 


137.1 


125.1 


112.6 


128.7 




Fc I 


5001.870 


3.881 


0.050 


102.5 




92.2 


109.0 




Fc I 


5014.951 


3.943 


-0.303 


102.9 




81.2 






Fc I 


5049.827 


2.279 


-1.355 




157.2 


134.8 


111.7 




Fe I 


5051.640 


0.915 


-2.764 




158.7 


144.0 


124.4 




Fe I 


5068.771 


2.940 


-1.041 


132.4 


111.5 


106.9 


85.8 




Fe I 


5074.753 


4.220 


-0.160 


85.2 




87.1 






Fe I 


5083.345 


0.958 


-2.842 




131.4 


122.3 


109.3 




Fe I 


5110.435 


4.260 


-3.758 








136.0 




Fe I 


5123.730 


1.011 


-3.058 








148.4 




Fel 


5127.368 


0.915 


-3.249 




113.9 








Fe I 


5216.283 


1.608 


-2.082 


142.5 


122.4 


140.3 


95.9 




Fel 


5225.534 


0.110 


-4.755 


114.2 










Fe I 


5232.952 


2.940 


-0.057 








141.9 




Fe I 


5269.550 


0.859 


-1.333 










156.0 


Fc I 


5281.798 


3.038 


-0.833 


142.6 




101.9 


75.1 




Fc 1 


5283.629 


3.241 


-0.524 




148.2 








Fe I 


5383.380 


4.312 


0.645 


126.4 


97.9 


103.7 


94.4 


34.0 


Fe I 


5393.176 


3.241 


-0.715 


117.1 


103.3 


91.0 


74.3 




Fe I 


5397.141 


0.915 


-1.982 










99.8 


Fe I 


5405.785 


0.990 


-1.852 










108.9 


Fe I 


5424.080 


4.320 


0.520 


149.5 


108.2 


111.9 


115.4 


39.2 
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Species 


A 


E.P. 


log gf 


EW(mA) 


EW(mA) 


EW(mA) 


EW(mA) 


EW(mA) 




(A) 


(cV) 




G108 


G322 


G315 


G351 


G219 


Fe I 


5429.706 


0.958 


-1.881 










122.5 


Fe I 


5434.534 


1.011 


-2.126 










98.7 


Fe I 


5446.924 


0.990 


-3.109 










103.7 


Fe I 


5497.526 


1.011 


-2.825 




143.7 




135.7 


66.6 


Fe I 


5501.477 


0.958 


-3.046 


147.7 


134.1 


122.2 


97.5 




Fc I 


5506.791 


0.990 


-2.789 




155.3 




136.0 


65.6 


Fe I 


5569.631 


3.417 


-0.500 


128.1 


104.9 


101.2 




30.3 


Fe I 


5572.851 


3.396 


-0.275 




141.8 


141.9 




38.9 


Fe I 


5576.099 


3.430 


-0.900 


103.6 


85.3 


81.3 






Fe I 


5586.771 


4.260 


-0.096 






145.2 




57.5 


Fe I 


5763.002 


4.209 


-0.450 


74.3 


73.5 


85.5 






Fe I 


6136.624 


2.453 


-1.410 










67.0 


Fe I 


6137.702 


2.588 


-1.346 




150.5 




105.5 


54.6 


Fe I 


6151.623 


2.180 


-3.330 


51.6 




47.2 






Fe I 


6173.341 


2.220 


-2.863 


71.7 


61.7 


60.9 


54.6 




Fe I 


6180.209 


2.730 


-2.628 






50.3 


26.0 




Fe I 


6187.995 


3.940 


-1.673 


40.7 




26.7 






Fc I 


6200.321 


2.610 


-2.386 


73.8 


54.5 


62.1 






Fc I 


6219.287 


2.200 


-2.428 


99.1 




81.9 


69.4 




Fc I 


6229.232 


2.830 


-2.821 


40.9 




29.9 






Fc I 


6230.736 


2.559 


-1.276 








125.2 


44.7 


Fe I 


6246.327 


3.600 


-0.796 


101.0 




91.9 


91.7 




Fe I 


6252.565 


2.404 


-1.767 


120.9 




104.3 


107.5 


44.5 


Fe I 


6254.253 


2.280 


-2.435 


114.4 




105.7 


82.6 




Fc I 


6265.141 


2.180 


-2.532 


96.8 




94.8 


72.4 




Fe I 


6270.231 


2.860 


-2.543 


46.0 




43.7 


22.1 




Fe I 


6297.799 


2.220 


-2.669 


101.4 










Fe I 


6335.337 


2.200 


-2.175 






108.9 


81.7 


31.0 


Fe I 


6336.830 


3.690 


-0.667 






89.9 


69.1 




Fe I 


6355.035 


2.840 


-2.328 






75.2 


51.8 




Fe I 


6393.612 


2.430 


-1.505 






129.3 


134.7 


47.3 


Fe I 


6411.658 


3.650 


-0.646 






94.4 


73.5 




Fe I 


6421.360 


2.280 


-1.979 








83.1 




Fe I 


6430.856 


2.180 


-1.954 






117.6 


86.3 




Fe I 


6481.878 


2.280 


-2.985 


74.3 










Fc I 


6494.994 


2.400 


-1.246 








151.1 


73.8 


Fe I 


6498.945 


0.960 


-4.675 


94.2 










Fe I 


6546.252 


2.750 


-1.536 


128.4 










Fc I 


6569.224 


4.730 


-0.380 


61.6 


52.7 








Fc I 


6593.874 


2.430 


-2.377 




85.1 


74.7 


64.5 




Fe I 


6677.997 


2.690 


-1.395 






123.1 




30.5 


Fe I 


6703.576 


2.760 


-3.059 


42.3 




39.9 






Fe I 


6710.323 


1.480 


-4.807 


44.3 










Fc I 


6750.164 


2.420 


-2.592 


89.7 


71.8 


71.1 


86.3 




re i 


OoUO.oOD 


9 7*^0 
Z. ( oU 


-Z.Ooo 












Fe I 


6839.835 


2.560 


-3.378 


28.8 










Fel 


6841.341 


4.610 


-0.733 






41.1 






Fe I 


7130.925 


4.300 


-0.708 


86.9 










Fe I 


7411.162 


4.280 


-0.287 




83.7 


72.6 






Fe I 


7445.758 


4.260 


0.053 


104.9 


92.2 


98.8 






Fe I 


7461.527 


2.560 


-3.507 


55.7 




45.3 






Fe I 


7491.652 


4.280 


-1.067 


70.1 


55.8 








Fe I 


7531.153 


4.370 


-0.557 






64.1 






Fe II 


4178.859 


2.583 


-2.489 




61.6 


51.1 






Fe II 


4178.859 


2.583 


-2.489 






78.5 






T?*:. TT 
re 11 


4Zoo. loy 


Z.DOO 


-l.yUU 






1 "^l 9 






Fc II 


4508.289 


2.856 


-2.318 


93.1 


73.2 


65.4 






Fe II 


4515.343 


2.844 


-2.422 


102.1 


84.3 


97.1 






Fc II 


4541.523 


2.856 


-3.030 






60.3 






T7o TT 
re 11 






1 son 
- 1 .oyu 


loU- 1 


±01-4l 


1 Af\ 9 


lOl -VJ 


fi^ 9 


Fe II 


4923.930 


2.891 


-1.307 


135.5 


118.8 


125.8 


121.8 


74.7 


Fe II 


5018.450 


2.891 


-1.292 




147.0 




136.7 


84.1 


Fe II 


5534.848 


3.245 


-2.790 


59.6 










Fe II 


6456.391 


3.903 


-2.075 






56.0 






Co I 


6770.970 


1.880 


-1.970 


49.5 










Col 


6814.961 


1.956 


-1.900 


45.8 










Co I 


6872.440 


2.010 


-1.850 


51.9 










Ni I 


6586.319 


1.951 


-2.810 


56.2 










Ni I 


6643.638 


1.676 


-2.300 


97.8 


97.7 


83.3 






Ni I 


6767.784 


1.826 


-2.170 


91.4 


79.3 


67.5 






Ni I 


7122.206 


3.542 


0.040 


98.6 










Ni I 


7393.609 


3.606 


-0.270 






66.8 
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TABLE 13 — Continued 



Species 


A 


E.P. 


log gf 


EW(mA) 


EW(mA) 


EW(mA) 


EW(mA) 


EW(mA) 




(A) 


(cV) 




G108 


G322 


G315 


G351 


G219 


Ni I 


7414.514 


1.986 


-2.570 


101.8 




63.3 






Ni I 


7422.286 


3.635 


-0.140 




63.3 


62.8 






Ni I 


7525.118 


3.635 


-0.520 






37.1 






Y II 


4883.690 


1.084 


0.070 


72.0 




73.5 






Ba II 


4554.036 


0.000 


0.163 










85.0 


Ba II 


4934.095 


0.000 


-0.157 










90.5 


Ba II 


5853.688 


0.604 


-1.010 






87.8 


86.5 




Ba II 


6141.727 


0.704 


-0.076 




142.4 




112.9 


67.5 


Ba II 


6496.908 


0.604 


-0.377 








141.8 


59.4 



Note. — 



Lines listed twice correspond to those measured in adjacent orders with overlapping wavelength coverage. 



